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STIMULATED  RAMAN  EFFECT 


INTRODUCTION  AND  SUMMARY 


This  report  is  concerned  with  phenomena  which  may  limit  the  transmission 
of  a  laser  beam  through  the  atmosphere.  We  have  considered  two  salient  types  of 
limitations  and  the  specific  contributing  processes.  The  first  limitation  cor¬ 
responds  to  a  self-defocus ing  instability  of  the  laser  beam  which  can  result  in 
an  undesirable  decrease  in  the  delivered  power  density.  The  instability  results 
from  photon  absorption  leading  to  atmospheric  heating  and  consequently  to  a  local 
decrease  in  the  refractive  index;  the  result  is  a  lens  effect  which  causes  the 
beam  to  diverge.  The  second  limitation  to  beam  transmission  is  imposed  by  the 
stimulated  Raman  effect.  This  process  involves  inelastic  photon  scattering,  and 
at  sufficiently  large  power  levels  results  in  the  amplification  of  the  beam  of 
scattered  photons  in  a  laser-like  fashion.  The  amplification  of  the  scattered 
photon  beam  can  produce  a  serious  depletion  of  the  laser  beam. 

Chapters  1,  2,  and  3  are  concerned  with  the  heating  mechanisms  that  lead 
to  the  self-defocusing  instability.  In  Chapter  1  the  far  wings  of  pressure- 
broadened  molecular  absorption  lines  are  studied.  The  motivation  is  the  need 
for  a  reliable  theory  to  compensate  for  the  absence  of  experimental  data.  The 
distance  from  the  line  center  beyond  which  the  usual  pressure-broadening  theory 
becomes  invalid  has  been  ascertained,  and  it  has  been  found  that  beyond  this 
point  the  line  profile  falls  exponentially.  This  conclusion  has  particular 
significance  with  regard  to  an  earlier  calculation  of  the  absorption  in  the 
window  at  1.0468  microns,  where  the  assumption  of  Lorentzian  line  shapes  led 
to  the  conclusion  that  the  absorption  coefficient  was  just  at  the  threshold 
for  the  self-defocusing  instability.  In  the  aforementioned  calculation  it  was 
assumed  that  the  absroption  is  due  only  to  the  molecular  resonance  lines.  The 
present  research  indicates  that  because  of  the  exponential  fall-off,  the  ab¬ 
sorption  due  to  the  molecular  resonance  lines  is  actually  substantially  smaller, 
and  consequently  may  be  ignored  with  respect  to  the  self-defocusing  instability 
at  1.0468  microns.  However,  it  has  been  pointed  out  that  atmospheric  absorption 
in  the  spectral  neighborhood  of  1.0468  microns  is  -actually  dominated,  not  by 
contributions  from  the  wings  of  molecular  resonance  lines,  but  by  a  continuous 
absorption  band  due  to  the  existence  of  ^2^2  comP^exe8,  Earlier  inferences 
from  experimental  data  indicated  that  the  absorption  coefficient  due  to  (^2^2 
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is  about  two  orders  of  magnitude  above  the  threshold  for  the  self-defocusing 
effect.  In  Chapter  2  we  report  the  results  of  a  detailed  study  of  the  photon 
absorption  mechanism  of  the  (02)2  complex  itself  and  the  relaxation  mechanisms 
by  which  the  absorbed  energy  is  subsequently  converted  into  heat.  It  is  con¬ 
cluded  that,  during  a  laser  pulse  of  about  one  millisecond  duration,  only  one- 
sixth  of  the  absorbed  photon  energy  is  converted  into  heat  so  that  the  estimate 
of  the  effective  absorption  coefficient  for  (^2 ^ 2  at  1*0468  microns  must  be 
reduced  by  a  factor  of  six.  Moreover,  the  study  suggests  that  the  desired  ad¬ 
ditional  reduction  in  the  net  absorption  coefficient  can  be  accomplished  by  a 
shift  to  slightly  higher  frequencies,  i.e.,  somewhere  in  the  region  1.030  to 
1.037  microns. 

Another  process  which  contributes  to  the  residual  absorption  in  what  is 
otherwise  considered  to  be  a  transmission  window  is  collision-induced  absorption. 
This  Is  a  phenomenon  in  which  absorption  takes  place  by  virtue  of  a  transient 
dipole  moment  which  exists  only  during  an  intermclecular  collision.  It  is 
characterized  spectrally  by  broad  continuous  bands  having  widths  of  the  order  of 
several  hundred  wave  numbers.  These  bands  are  superposed  upon  internal  transition 
lines  which  may  even  be  +n  bidden  for  the  isolated  molecule. 

The  energy  of  th.  absorbed  photon  is  partitioned  into  an  internal  ex¬ 
citation  of  one  of  the  molecules  and  kinetic  energy  of  the  colliding  pair.  The 
energy  taken  up  by  the  translational  component  corresponds  to  an  instantaneous 
heating,  which  contributes  to  the  beam  instability  regardless  of  the  laser  pulse 
duration.  The  translational  band  contributes  to  the  magnitude  of  the  absorption 
in  a  valley  between  resonance  lines,  and  may  dominate  over  the  far  wing  con¬ 
tributions  of  resonance  lines. 

In  Chapter  3  we  consider  the  photon  absorption  which  occurs  during  a  col¬ 
lision  between  two  rare  gas  atoms.  The  analysis  is  somewhat  simplified  in  this 
case  because  of  the  absence  of  rotational  and  vibrational  modes;  the  absorbed 
energy  goes  entirely  into  translational  energy  of  the  atoms.  Theoretical  work 
in  this  field  has  been  meager  and  has  been  principally  confined  to  the  calcu¬ 
lation  of  frequency  integrated  absorption  coefficients.  Relatively  little  has 
been  done  regarding  the  frequency  dependence  of  the  collision  induced  absorption 
coefficient  and  even  this  work  is  quite  recent.  The  studies  which  have  been 
made  were  not  concerned  with  the  far  wings  and  cannot  be  extrapolated  into  the 
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contains  a  discussion  of 


far  wing  regions  with  confidence;  the  present  chapter 
these  calculations. 

We  have  performed  a  quantum  mechanical  analysis  of  this  problem  vhich  Is 
capable  of  yielding  the  absorption  coefficient  as  a  function  of  frequency  even 
in  far  lings  of  the  absorption  band,  where  no  experimental  determine  on^as 
as  yet  been  possible.  It  Is  these  unmeasurable,  far  wing  absorptions  which 
tend  into  the  atmospheric  windows,  and  provide  Instantaneous  heating.  The,  mu. 
be  determined  theoretically  in  order  to  evaluate  the  self-defocusing  heating  in¬ 
stability.  The  mathematical  analysis  leading  to  the  expression  for  the  absorp  on 
efficient  has  been  completed,  and  only  a  number  of  independent  numeric,  inte¬ 
grations  remain  to  be  performed.  Preparations  are  currently  in  progress  for 
evaluation  by  digital  computer. 

in  Chapter  4  laser  beam  depletion  by  stimulated  Raman  scattering  Is  in¬ 
vestigated  by  solving  rate  equation,  describing  this  process.  The  “ 

effects  whil  are  analysed  in  great  detail  are  (a,  the  interaction  .  «-, ^ 

Raman  lines,  (b)  multiple  order  Stokes  wave  generation,  (c)  the  amp  ca 

,  ..  ...  anA  /d\  the  effect  of  pressure  dependence  upon  the 
Raman  backscattered  light,  and  (d)  the  erre 

propagation  of  a  beam  vertically  through  the  atmosphere.  In  each 

the  intensities  of  both  the  Raman  radiation  and  primary  laser  beam  are  o  an 

as  a  function  of  distance  from  the  laser. 

It  is  found  that  the  maximum  intensity  of  laser  light  which  can  be 
mi t ted  vertically  through  the  atmosphere  without  frequency  change  due  to  ma 
scattering  is  critically  sensitive  to  the  laser  line  width.  e  max  m 
tensity  is  calculated  for  representative  values  of  the  laser  ^ 

other  pertinent  parameters.  Another  conclusion  is  that  the 
scattered  beam  has  a  significant  deleterious  effect  if  the  ° 

backscattered  cross  section  is  greater  than  ninety-two  percent  of  the  forward 
cross  section,  and  the  laser  pulse  duration  is  greater  than  5  x  10  secon  s. 

I£  either  of  these  conditions  fails  to  prevail,  then  Raman  backscattering 

inconsequential  as  a  loss  mechanism. 

An  analysis  is  reported  which  describes  the  propagation  of  a  pulse  of 
radiation  moving  through  and  antiparallel  to  the  laser  beam.  The  frequency 
the  backward  propagating  pulse  is  that  which  corresponds  to  the  Stokes  shifted 
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frequency.  It  is  found  that  this  pulse  may  be  amplified  at  the  expense  of  the 
laser  beam  to  a  peak  intensity  many  times  greater  than  that  of  the  laser  beam 

itself,  and  that  in  fact  most  of  the  energy  in  the  laser  beam  can  be  concentrated 
into  this  backward  moving  pulse. 


( 
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Chapter  1 


FAR  WING  PRESSURE  BROADENING 


In  this  chapter,  some  model-type  studies  on  the  far  wing  pressure  broadening 
of  molecular  vibrational  lines  will  be  reported.  The  most  elementary  system  of 
interest  is  that  of  a  diatomic  heteropolar  molecule  in  collision  with  a  monatomic 
perturber.  For  the  sake  of  definiteness,  the  optical  transition  will  be  assumed 
to  connect  the  vibrational  state  with  the  first  excited  (vibrational)  state  of 
the  diatomic  molecule  (0  ♦+  1) . 


The  physical  picture  motivating  the  treatment  derives  from  earlier  work'L 
by  one  of  the  present  authors;  it  may  be  described  briefly  as  follows.  When  the 
optical  frequency,  w,  is  sufficiently  close  to  the  unperturbed  transition  fre¬ 
quency,  —  specifically  when  (w- a)0)Tc<  with  ic  an  appropriately  defined 
collision  time,  one  may  expect  the  impact  theory  to  be  valid;  the  line  shape  will 
therefore  correspond  to  the  wing  of  a  Lorentz-type  line  shape  7^/(0)- wq)  ^ .  If, 
on  the  other  hand,  the  frequency  shift  (uj— )  is  so  large  that  (w-u)o)tc>  1,  the 
spectral  profile  will  follow  the  predictions  of  the  so-called  statistical  theory. 
In  this  theory,  one  focuses  attention  on  the  collision-induced  perturbation  of 
the  instantaneous  transition  frequency.  Two  possibilities  then  arise.  If  during 
a  collision  event,  the  perturbed  transition  frequency,  u)q+  momentarily  co¬ 
incides  with  the  optical  frequency,  oj,  the  spectral  intensity,  I(cu),  will  be 
proportional  to  the  occurrence  probability  of  the  perturbed  transition  frequency 
—  i.e.,  the  probability  that  ujo+  u)p  lie  within  a  unit  differential  frequency 
range  about  the  spectral  frequency,  a>.  If,  on  the  other  hand,  there  are  no 
collisions  which  are  able  to  provide  a  coincidence  of  ojo+  with  w,  one  may 
expect  the  spectral  intensity,  I(oj),  to  decrease  exponentially^  with  increasing 
frequency  shift,  w-oj  . 


The  particular  problem  upon  which  the  present  studies  are  focused  is  the 

absorption  on  the  violet  side  of  a  given  vibrational  line,  at  distances  of  the 
3  -1 

order  of  600-1000  cm  from  the  unperturbed  frequency.  (The  corresponding 

14  14  -1 

shift  in  units  of  circular  frequency  is  u>-o)o~  1.2x10  —  2x10  sec  .)  From 

the  above  discussion  it  then  becomes  clear  that  the  absorption  at  such  extreme 
distances  from  the  unperturbed  line  will  be  exponentially  small  unless  (a)  the 
collision  time  is  <  10  ^sec  ^  —  with  atomic  velocities  **'10^cm/sec,  such  3 

_9 

limitation  could  be  met  only  with  very  short  range  interactions  (~10  cm) 
and/or  (b)  the  coll  is  ion- induced  shift  in  the  momentary  vibrational 
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frequency  is  positive  and  Z  600-1000  cm  Whether  or  not  such  extreme  conditions 
can  be  met,  one  may  infer  that  they  will  be  approached  most  closely  by  the  standard 

A 

short-range  repulsive  interatomic  interactions  generally  taken  to  have  the  form 

V(r )  «  V~ar  ’  (1.1) 

rather  than  the  longer  range,  multipole-type  interactions  of  the  form,  C/rn 
(with  n  a  number  varying  roughly  from  3  6).  It  is  therefore  proposed  to  consider 

the  basic  interaction  responsible  for  the  far-wing  absorption  to  be  a  sum  of  two 
terms,  each  of  the  form  (1.1)  —  describing  the  interaction  between  the  colliding 
atom  and  each  of  the  two  atomic  constituents  of  the  absorbing  molecule.  Moreover, 
in  order  to  avoid  complications  which  —  in  the  opinion  of  the  present  authors  — 
do  not  appear  to  be  essential  for  the  physics  of  the  problem,  the  following 
further  "model-type"  simplifications  are  introduced: 

(1)  the  rotational  motion  of  the  molecule  is  ignored,  and 

(2)  the  collision  is  assumed  to  be  collinear. 

Concerning  the  first  simplification,  it  may  be  remarked  that,  on  the  scale 
of  the  above  mentioned  frequency  shift  (600  -  1000  cm  ^) ,  rotational  frequencies 
are  quite  small,  so  that  (a)  there  is  certainly  no  question  of  the  existence  of 
any  rotational  component  of  the  line  in  the  far-wing  spectral  region  of  interest 
and  (b)  ,  in  considering  the  partition  of  the  excess  energy  into  translational  and 
rotational  motion,  it  would  seem  quite  legitimate  to  ignore  the  discrete  structure 
of  rotational  energy  levels.**  Once  this  is  done,  rotation  and  translation  may 
both  be  considered  as  comparable  types  of  essentially  unquantized  motion. ^  For 
a  detailed  quantitative  calculation,  both,  of  course,  must  be  taken  into  account. 
For  a  semiquantitative  study  which,  like  the  present  one,  is  aimed  at  a  first 
look  at  the  gross  features  of  the  far-wiq,g  broadening  problem,  it  is  felt  that 
a  model-type  treatment  which  takes  explicit  account  of  energy  conversion  into 
one  type  of  unquantized  motion,  namely  translation,  while  ignoring  the  rotational 
component,  is  basically  adequate. 

The  assumption  of  collinear  collisions  is  introduced  because  it  is  sus¬ 
ceptible  to  a  relatively  straight-forward  analysis.  It  is  hoped  that  the  insight 
gained  in  studying  this  special  case  will  permit  at  least  semiquantitative  ex¬ 
tensions  to  the  more  general  case. 
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It  may  finally  be  stated  that  both  of  the  above  discussed  simplifications 
are  contained  in  essentially  all  treatments  of  the  related  problem  of  collision- 
induced  vibrational  deactivation  of  diatomic  molecules. 


With  these  preliminary  considerations  out  of  the  way,  let  us  proceed  to 
the  formulation  of  the  model  Hamiltonian  of  the  problem.  In  the  absence  of  the 
external  (optical)  electromagnetic  field  it  takes  the  form 


p2  p2 

h  -  Hr  *  i  “v“„2  x  +  IT  +  v(p-x> 


(1.2) 


In  this  expression  X  £  r  -  rg  (with  r  the  internuclear  separation  and  rg  its 

equilibrium  value)  is  the  vibrational  coordinate,  px  *  J  the  associated 
canonical  momentum,  -  *  the  reduced  mass  for  vibrational  motion  (letters 

A  and  B  will  be  used  lo  designate  quantities  associated  with  the  constituent  atoms 
of  the  diatomic  molecule,  whereas  the  letter  C  will  refer  to  the  colliding  mon¬ 
atomic  species) ,  R  denotes  the  distance  between  the  center  of  gravity  of  the 
molecule  and  the  collidine  atom  (C) ,  Pg  .  *  &  1-  the  associated  canonical 

momentum  (with  M  .  the  reduced  mass  for  the  relative  motion  of  atom  C 

with  respect  to  the  cintlr  of  gravity  of  the  molecule).  The  first  two  terms  in 
(1.2)  are  the  kinetic  and  potential  energy  for  vihrational  motion  of  the  molecule 
(„  .  Vibrational  frequency).  The  third  term  is  the  kinetic  energy  for  the 

relative  motion  of  the  colliding  partners.  Finally.  V(R.X)  Is  the  interaction 
between  the  colliding  atom  C  and  the  molecule.  It  is  here  assumed  to  have  the 

form 


V(R,X)  =  Vq  exp{-a[R-*r]}  *  Vq  exp{-a[R-A(re+  X) ] } 


(1.3) 


where 

»  =  -A-  .  «-4> 

■  «a+mb 

Expression  (1.3)  follows  from  (1.1)  and  the  fact  that,  in  a  collinear  collision, 

the  interaction  of  colliding  atom  C  with  only  ore  atom  of  the  molecule  (say,  B) 

need  be  considered  (cf  numerical  estimate  of  a  given  in  Ref.  4  and  the  tables  of 

values  of  r  in  Herzberg's  book  ). 
e 
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The  presence  of  the  optical  electromagnetic  field  gives  rise  to  an  ad¬ 
ditional  term  in  the  Hamiltonian.  It  has  the  form 

H  -  -y(X)  £jfia)t  +C-C*  (1.5) 

em 

in  which  £  is  the  amplitude  of  the  electromagnetic  field  of  the  incident  light 
and  y(X)  is  that  part  of  the  molecular  dipole  moment  which  depends  upon  the 
vibrational  coordinate,  X;  as  is  usual,  it  will  he  assumed  that  this  dependence 
is  linear,  so  that  the  optical  selection  rule  for  the  change  in  vibrational 
quantum  number  Av  ■  1  1. 

In  the  absence  of  the  electromagnetic  perturbation,  represented  by 
(1.5),  there  exist  stationary  solutions,  ij^(X,R)  of  the  time-independent 
Schroedinger  equation 

H  ^(X.R)  =  E1*1(X,R)  ,  (1.6) 

which  are  asymptotically  (R  •+  ®)  of  the  form 

<J>.  (X,R)  -  <f>  (X)  sin(k,R+<5),  (1.7) 

l  o  l 

where  <f>Q(X)  is  the  ground  state  Harmonic  oscillator  wave  function  and 

k^  =  /2ME^/tl^  is  the  incident  wave  vector  of  relative  motion.  In  what  follows, 

iHX,R)  will  be  approximated  by  a  product  function 


<J»  (X,R)  •*  MX)  MR), 
i  o  o 


where  ^(R)  is  that  solution  of  the  equation 


M  (R) 
1  o 


^R 

2M+Voo(R)|  VR)’ 


0 


v  (R) 

oo 


$  (X)  V(R, X)  dx  , 
o 


which  obeys  (1.7). 


(1.8) 


(19) 

(1.10) 
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In  the  presence  of  the  electromagnetic  field,  the  wave  function  of  the 
system  will  be  augmented  by  a  term,  ^(X.R),  which  may  be  computed  from  first- 
order  time-dependent  perturbation  theory.  The  relevant  equation  is10 

3*1 

9t“  “  H*1  "  Hem  *i(x»R>  exP{  (it/1i)(E1+tu1Jo/2)j,  (1.11) 


W±th  Hem  glv®"  by  The  steady-state  solutions  of  (1.11)  will  consist  of 

two  terms,  *  and  ^  with  time-dependences  of  the  form 

ejtp{-(it/il)(E1+‘naio/2  +  1ki)}  and  |-(it/ft)  (E±+  1iu>o/2^H(o) }  respectively;  of  these, 
only  the  former  (describing  the  absorption  of  a  photon  of  energy  •«(*)  is  of  interest). 
The  equation  determining  this  term  is 


■ftw 

+  2°  +  ^ 


-U(X)Ct|,i(X,R). 


(1.12) 


At  this  point,  a  considerable  simplification  can  be  achieved  by  resorting 
to  the  (already  employed)  distorted-wave  approximation.  This  may  be  do-«.  as 
follows.  Substituting  the  expansion 


*<+  "  2  Xm+)*m(X)  ’  (1.13) 

m 

f 

[where  the  *  (X)  are  the  eigenstates  of  the  vibrational  part,  H  ■  — +  —M  “  ^*2 
“  I  1\)  vib  2\i  2  vo  * 

with  eigenvalues  ‘ft^Jn  +  ^Jj  multiplying  on  the  left  by  *n(X),  and  integrating 

with  respect  to  X,  one  has  [upon  using  (1.2)  and  (1.8)] 


Ei+  2M  ' 


V  (R) 

nn 


)xn+>  (R)  +  I  V. 
’  myn 


nm  (R)X,1+)(R) 

nm  m 


~Zv  *  (R) , 
no  o 


(1.14) 


where 


nm 


V  (R) 

r,m 


[  *  (X)p(X)*  (X)  dX  , 

;  n  n 

|  *n(X)V(R,X)*m(X)  dX  . 


(1.15) 

(1.16) 
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The  distorted-wave  approximation  consists  in  ignoring  the  sum  on  the  l.h.s.  of 

(1.14),  (i.e.,  terms  containing  the  coefficients  V  ,  with  mj*n) .  One  then  has, 

nm 

in  place  of  (1.14) 


E .  +  1i  (to-u)  ) 
1  o 


Pr 

2M  -  V11(R) 


>x<+) 


(R) 


■*u10*o 


(R) 


(1.17) 


since,  with  p(X)  assumed  to  be  linear  in  X,  p  differs  from  zero  only  for  n  =  1. 

no 


It  is  now  convenient  to  write 


in  the  form 


. . +  x  (E)  . 
H(w-M  )  XSCW  * 
o 


(1.18) 


upon  substituting  (1.18)  into  (1.17)  and  making  use  of  (1.9),  one  obtains  for 

\ l>  (R)  the  equation 
sc 


E1+  H(U-Uo)  -  ^  - 


vu(R) 


h(o)-u)  ) 
o 


£uio^o(R)- 


(1.19) 


Before  proceeding  further,  some  comments  on  the  physical  significance  of 
the  two  terms  on  the  r.h.s.  of  (1.18)  are  in  order.  Briefly,  it  is  claimed  that 
the  second  term,  x  (R) ,  represents  the  amplitude  of  the  process  under  study  here 
—  namely,  the  energy-conserving  process  whereby  a  quantum  of  energy,  Ito,  is 
absorbed,  with  a  final  translational  energy  (of  relative  motion),  Ej+1l(u>-u>  ) 
(as  shown  by  the  form  of  the  l.h.s.  of  (1.19)).  By  way  of  contrast,  the  first 
term  on  the  r.h.s.  of  (1.18)  represents  the  virtual  absorption  of  a  quantum 
(in  a  0  -►  1  vibrational  transition);  it  particular,  from  (1.8)  one  sees  that, 
since  the  R  dependence  is  contained  entirely  in  the  factor  ^o(R),  the  kinetic 
energy  associated  with  this  state  is  the  initial  E^  [rather  than  E^  ■  E^+  (cu— u)q )  ] . 
The  virtual  character  of  the  transition  described  by  the  first  term  of  (1.18)  is 
thus  manifest;  it  will  therefore  be  discarded  without  further  ado. 

From  a  "technical"  standpoint,  it  may  be  remarked  that  Eq.  (1.19)  for 
Xgc(R)  is  rather  more  suitable  for  scattering  studies  than  the  corresponding 

i 


it 


Eq.  (1.17)  for  (R)  •  The  reason  Is,  basically,  that,  in  contrast  to  (1.17), 
the  driving  term  of  (1.18)  contains  the  factor  ^V^(R)  -  Vqo(R)J  ,  which  differs 
from  zero  only  in  the  immediate  neighborhood  of  the  collision  region  (R  ** 


interatomic  distance).  Outside  this  region,  x  (R)  satisfies  the  homogeneous 


wave  equation;  the  usual  boundary  condition  of  scattering  theory  —  namely,  that 


as  R  4  ",  X  (R)  takes  the  form  of  an  outgoing  wave(“i  f  )  is  then  straight- 

SC 


forwardly  applied,  as  will  now  be  done. 


One  multiplies  both  sides  of  (1.19)  by  that  solution,  ij^(R),  of  the  homo¬ 
geneous  equation: 


PR  1 

Ef  -  2iT  V11(R)JV 


(1.20) 


which  has  the  asymptotic  form 


^(R)  -■--■"-♦sin 


[k£R«] 


(1.21) 


and  integrates  from  zero  to  R.  Performing  the  approximate  integrations  by  parts 

(  2  2  2  2\ 
with  PR  =  -ft  3  /3R  ) 


i(^-)  |VR’>[V11<R,>  -  Voo(R'>]  VR’>dR' 

n  *  n 


2f  d*.(R) 

likc(R)  -n~ 

-  *i<R)dr  v<R]  • 


(1.22) 


For  sufficient  large  R,  the  upper  integration-limit  on  the  l.h.s.  of  (1.21)  may 
be  replaced  by  infinity;  at  the  same  time,  one  may  introduce  (1.21)  and  the 


asymptotic  expression 


R  -»•  »  .  .ikcR 


J$r. 


imi 


tCr-' 


into  (1.22).  This  procedure  yields  sn  explicit  expression  for  the  scattering 


amplitude,  A  ,  viz: 
sc 


c^io  p 

•h(0)-0)o)  [V11(R)_V00(R)]  dR- 


(1.24) 


The  outgoing  particle  flux  is  given  by  the  relation 


^ut  "  vflAsc 


o  flic- 

I  ■  irl*scl2  • 


Introducing  (1.24)  Into  (1.25)  end  dividing  by  the  Incident  flux, 


(1.25) 


rin  "  V4  "  * 


(1.26) 


one  obtains  a  quantity 


_  4kf  4H2 


_f  4jr_  «•  "ip  r  ,  ,  ,2 

ki  1,V„)R  J  *l(R,*o(R,[Vll(,»-''oo(R)]d»  ■  (1.27) 


o'  1  O 


vhich  nay  be  regarded  as  the  probability  that,  in  a  collision  of  the  speciel 
type  under  consideretion  -  nanely,  a  unilinear  collision  -  a  (0  -  1)  vibrational 
transition  occurs  (with,  of  course,  sinultsneou.  ebsorptlon  of  a  quantum,  -hn, 
of  electromagnetic  energy). 

In  order  to  proceed  further,  it  is  desirable  to  write  down  explicit  ex¬ 
pressions  for  Vn(R)  and  V  (R) .  From  (1.3),  one  has 


V  (R)  .  V  e_aR  . 

oo  i  * 


(1.28a) 


V„(R)  -  V,e”aR 


fe  ’ 


(1.28b) 
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where 


14 


becomes 


aft  j  ist  + 

d52  sd5 


ijik2 

2  2 
•ot  5 


this  equation  is  to  be  supplemented  by  the  boundary  conditions 


(1.33) 


S-+-  CD 

- ►  0  , 


(1.34a) 


♦ 


->  sin 


2k 


2b 


—  lo«  a?  + 


(1.34b) 


The  required  solution  is 


♦  - 


J(2k/air)  sinh  2irk/aJ 


1/2 


K2ik/a(5)  ’ 


wheru 


K  (?) 
v 


sin  irv 


(1.35) 


(1.36) 


is  (in  standard  notation)  the  Bessel  function  of  imaginary  argument  which  obeys 
condition  (1.34a).  To  see  that  (1.34b)  is  satisfied  one  notes  that,  from  (1.36) 
and  Eq.  (2),  section  3.7  (p . 77)  of  Watson's  treatise  on  Bessel  functions, 


K2ik/a(5) 


2i 


.  .  2trk 
sinh  1 
a 


(2/Q 

r(i-2 


_  (2/e) 


-2ik/a 


2ik/a 


2ik7cO'  r(l+2ik/o) 


■  '(2k/a)  Blnh  »1"[(2k/“>  1"*<2^)+S']  •  (1'37) 

the  equality  following  from  the  standard  formula 

|r(l+iy)|2  -  iry/sinh  iry  » 
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[with  6'  =  arg  r(l+iy)].  From  (1.37),  it  is  clear  that  (1.35)  satisfies  (1.34b)* 


Eq.  (1.35)  is  now  utilized  to  provide  explicit  expressions 
*0(R>  and  ^(R)  in  Eq.  (1.27);  the  result  is 


10 


£ 


,2 


10 


•ft(o)-o)  ) 
o 


(W~  o  27rkf  2Trke  o 

H }'  v  (*/"2>  sinh  -£-*■  sinh  -~|^|2  » 


where 


K2ik4J/a(z)zK2ik1/a(6z)dz  » 


o  f 


with  (cf  (1.30c)  and  Ref.  11) 


.2  =  bi  Vt  i+a2*^/Hwvu0 


b2  Vf  l+3a2X^ti/^ 


y  u> 

v  o 


It  is  indeed  fortunate  that  the  integral  on  the  r.h.s. 
explicitly  evaluated.  Specifically,  one  has  [cf  Tables  of  In 
Bateman  Manuscript  Project,  Vol.  I,  p.  334,  No.  (47)]. 


of  C 

egra! 


6  |  Kjj(z)Kv(6z)z{:z 

*  O 


1 

2 


8*1 

r(2) 


X 


F 

2*1 


U+V 

2 


+  1, 


2; 


where  2Fl^a,b,c;z^  is  the  hyper geometric  function  which,  for  |z|  < 
by  the  series 


2*2.  ^®»bic;z) 


r  (c)  p  r(a+n)r(b+n) 

r(a)r(b)  L  r(c+n)  n'. 
n-o 


n 

z  . 
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x 


for  both 

(1.38) 

(1*39) 

(1*40) 

.39)  has  been 
Transforms. 

(1.41) 

| ,  is  defined 

(1.42) 


It  possesses  the  Integral  representation  (cf  Handbook  of  Mathematical 
Functions .  Applied  Mathematics,  Series  55,  U.S.  Department  of  Commerce,  National 
Bureau  of  Standards,  p.SSRj  Frj.  15.3.1) 


2^i 

,c;z)  - 

r(c)  f1  tb_1(i-t)c"b'1 

(1.43) 

r(b)r(c-b)  Jo  (1.tl). 

valid  when 

Rc  > 

Rb  >  0  . 

(1.44) 

In  the  case  at 

hand 

a  ■ 

2HL+  1 

2  A 

» 

(1.45a) 

b  - 

9 

(1.45b) 

C  “ 

2  , 

(l.45c) 

z  » 

i-e2  , 

(1.45d) 

and  one  notes  that  [cf  (1.39)  and  (1.41)]  with  u  and  v  imaginary,  (1.44)  is 
obeyed.  Substituting  (1.43)  and  (1.45)  into  (1.41),  one  then  has 


ej  Ky(z)Kv(Bz)zdz  -  j  Bv+1r 


i 


i_ 

2 


X 


dt 

[l-t(l-B^ 


(u+v)/2  +  1 


9 


which,  by  virtue  of  the  relationship 

r(l+z)T(l-z)  -  zff/sin  irz  , 
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may  be  written  as 


8 


Ky(z)Kv(gz)zdz 


sin  (p+v)ir/2 


1  ftv+l 

2  6 


r 


dt 


[l-t(l-S2j 


<U+v)/2  +  l’ 


(1.46) 


It  is  expedient  to  transform  the  r.h.s.  of  (1.46)  via  the  introduction  of 
a  jiew  integration  variable,  x  -  1/t;  the  result  is  [cf  (1.39)] 

.  1  ^(u-fvH/z  r  d»,<|1-K»2f,-n<>'-v>/2 

2  sin  ir(y+v)/2  J  .(p4v')/2  +  1  *  (1-47) 

x[x-(l-62)j 

where 


li  =  2ikf/o  , 


(1.48a) 


and 


v  =  2ik^/a  . 


(1.48b) 


The  stage  has  now  been  reached  where  the  various  limiting  cases  outlined 
in  the  introduction  may  be  studied.!  Let  us  first  consider  the  domain  of  the 
statistical  theory.  From  the  previous  work1  (cited  in  the  introduction)  it  may 
be  inferred  that  this  domain  corresponds  to  those  values  of  the  parameters  p,  v, 
aad  0  for  which  (1.47)  may  be  evaluated  by  the  method  of  steepest  descents.  Let 
us  therefore  carry  out  such  an  evaluation;  besides  yielding  concrete  results  for 
the  statistical  limit,  the  treatment  will  serve  to  delineate  the  range  of  values 
of  u,  v,  and  0,  for  which  this  limit  constitutes  a  good  approximation. 
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descent  evaluation  of  (1.47),  it  is 


For  the  purpose  of  the  steepest 
convenient  to  write  it  in  the  form 

•u  1  Rv41(u+v)tt/2  [“  d  x  eF(x)_ 
^  “  2  sin  ir(y+v)/2  x(x-c) 

where 

c  =  1  -  82  , 


(1.49) 


(1.50) 


and 

F(x)  ■  ^2^'  1°8  x  +  1°£>  (x“^  "  ^ 2  ^x"c^ 

=  [log  x  +  Y  log  (x-1)  -  log  (x-c)J  , 

with 


y+v 


(1.52) 


Under  the  conditions 


|Ji±£|  »  l  , 

(1.53a) 

I  j  1 

p—i  »  i . 

(1.53b) 

I  j  1 

; 

the  factor  eF(x)  is  a  rapidly  varying  function  of  x, 

except  in  the  vicinity  of 

the  generally  complex  saddle  point,  zg,  defined  by  the  relation 


(1.54) 


(1.54)  constitutes  a  quadratic  equation  in  the  unknown  saddle  point,  zg. 


X 
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The  two  solutions  are 


z 

s 


+ 


(1+1/Y)2  -  c/y 


(1.55) 


The  task  of  demonstrating  that  the  original  integration  contour  of  (1.49) (namely, 
along  the  real  axis  from  1  to  “)  may  be  deformed  into  one  whose  principal  con¬ 
tributions  occur  at  the  saddle  points  is  relegated  to  a  future  report.  Preliminary 
studies  show  that  (a)  if  the  saddle  points  [given  by  (1.55)]  are  complex,  the 
major  contribution  to  (1.49)  occurs  along  a  steepest-descent  contour  in  the  im¬ 
mediate  vicinity  of  that  saddle  point  which  is  located  in  the  upper  half  of  the 
complex  plane,  and  (b)  if  the  saddle  points  are  real  and  greater  than  unity 
(this  latter  property  holding  when  c  <  1,  as  will  hereinafter  be  assumed  to  be 
the  case),  contributions  to  (1.49)  occur  in  the  neighborhoods  of  both  saddle 
points.  In  both  cases,  the  evaluation  of  these  contributions  requires  computing 
F"(zs).  For  algebraic  convenience,  it  is  desirable  to  write  (1.54)  in  the  form 


F'(z  ) 
s 


where 


M±L 

2z 

s 


g(l/z8) 


0  , 


g(y) 


l 


(1.56) 


(1.57) 


One  then  has  (upon  using  (1.56)  and  introducing  the  notation  y  ■  1/z  ) 

8  8 


F"  (z  ) 
s 


p+v 

2 


^7 


c 

(l-cy8)2 


The  solution  of 


(1.56)  in  terms  of  y 

Js 


reads 


where 


-  i±X  + 
2  - 


(1.58) 


(1.59a) 


(1.59b) 
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Upon  introducing  (1.59a),  (1.58),  and  eliminating  c  by  the  use  of  (1  59b),  one 
has  (after  some  algebraic  manipulation) 


At  this  point,  it  may  parenthetically 
the  factor  z^z_c^  ,  present  in  the  integrand 
brought  into  the  form 


be  noted  that,  by  similar  manipulation, 

of  (1.49),  may  (for  z  ■  z  )  be 

8 


(1.61) 


The  insertion  of  (1.60)  and  (1.61)  into  the  expression  for  the  saddle  point 
integrals  yield  (as  will  be  seen  immediately  below)  numerically  equal  con¬ 
tributions  when  the  saddle  points  are  on  the  real  axis,  a  feature  which  will 
receive  physical  interpretation. 


In  the  case  of  real  saddle  points  (A  real)  F"(zs)  [as  given  by  (1.60)] 
is  manifestly  imaginary  (y,  and  v  are  imaginary);  hence  the  contour  of  steepest 
descents  is  inclined  at  an  angle  of  ir/4  or  3ir/4  with  respect  to  the  real  axis 

[depending  on  the  sign  of  F"(z)].  The  saddle  point  integrals  are  then  of  the 
form  [cf . (1.49)  ] . 


V  -  1  exP(*l*/4)  exp(F(zs)) 

s  2  sin(y+v)ir/2  ' 


Z8(ZS-C) 


+» 


exp(  |F"(z  )  |r>/2)  dn 


1  BV+1(h+v)7t/2  e*P(*l*/4)  exp(F(zfl)) 

2  sin(p+v)7r/2  z  (z  -c) 

x  s 


2?r 


I F" (zs) | 

which,  with  the  use  of  (1.60)  and  (1.61),  may  be  written  as 


1/2 


"  2  c^sin(u+v)Tr/2  exP(±llr/4)  ^P^)) 


2-nv 


(y+v)  A 


1/2 


(1.62) 
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Bearing  in  mind  that  (a),  F(zfi)  is  Imaginary  for  real  saddle  points  and 

(b)  in  an  eventual  energy  integration,  the  relative  phase  of  F(z  )  at  the  two 

s 

saddle  points  will  undergo  many  oscillations,  let  us  simply  ignore  the  inter¬ 
ference  term  which  results  from  the  substitution  of 


x-T*.. 

s-r,2 

into  (1.38).  One  then  has 


(1.63) 


r 


10 


10 


fi(u-u)  ) 
o 


V  V 
f  i 


(1.64) 


where,  with  use  of  (1.50),  (1.48a,b),  (1.40),  and  (1.59b),  one  obtains  (after 
some  algebraic  manipulation) 


1 

2 


8  |  Iu+vItt^/4 _ 

■1-82'  sinh2(  |y+v  |tt/2) 


1 

A 


V  V 

i  Vf 


.2  .  2 

kf~ki 


2  (Vf'Vi)2 


tt3/4 


*lnh2[a(kf+kl)] 


where 


2 

ic 


k  -k 
.2  Kf  Ki 

Kf  ’  i-yvf  ' 


(1.65) 


(1.66) 


A  brief  digression  on  the  physical  significance  of  the  quantity,  k  is  now 
in  order.  From  Eqs.  (1.30a,b,c)  and  (1.28a,b),  one  notes  that 


and 


<i  (R) 


2 

i 


<fW  = 


.2,2  -aR 
kf-bf  e 


,1/2 


(J. 67a) 


(1.67b) 
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1 


are  "local"  wave  vectors  for  translational  motion  along  the  potential  curves 
V  (R)  and  Vn(R)  respectively,  (with,  of  course,  total  energies  Ej^  and  Ef^ 
respectively).  The  condition  that  these  two  local  wave  vectors  shall  coincide 

at  some  common  value,  k,  i.e., 


kl(R) 


k^(R)  =  k  t 


(1.63) 


is  readily  shown  to  be  equivalent  to  (1.66).  It  is  thus  seen  that  the  principal 
contributions  to  (1.65)  [i.e.,  to  the  matrix  element  for  the  transition,  as 
given  by  (1.27)]  obey  a  modified  Franck-Condon  Principle,  in  which  the  instantaneous 
position  and  velocity  of  relative  translational  motion  are  conserved. 

The  above  remarks  lead  to  a  simple  physical  interpretation  for  the  occur¬ 
rence  of  the  factor,  k,  in  the  denominator  of  (1.65).  Bearing  in  mind  that 


Mv 


in 


(1.69) 


(where  v  is  the  common  instantaneous  velocity  of  relative  motion  in  the  initial 
and  fina^states),  one  sees  that  the  1/k  dependence  is  simply  an  expression  of 
the  usual  1/v  law  for  the  occurrence-probability  distribution.  This  inter¬ 
pretation  will  be  developed  more  explicitly  immediately  below. 

Proceeding  with  the  calculation  one  has,  upon  inserting  (1.65)  into  (1.64), 
and  utilizing  the  conservation-of-energy  relation 


•ft(u)-w  )  *  , 


(1.70) 


together  with  (1.69)  one  has 


e2  2 

c  u10 


2* 


10 


2  v.a(u)-u)  ) 

in  o 


|l-exp  (-4itk^/ a)j  fl-exp  (-4itk^/ a) 
|l-exp(-2if  (k^+k^)  /a)J  ^ 


(1.71) 
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In  discussing  the  relationship  of  (1.71)  to  the  prediction  of  the  sta¬ 
tistical  theory,  let  us  first  note  that,  by  virtue  of  conditions  (1.53a,b) 
(which  define  the  domain  of  validity  of  (1.71)],  the  curly  bracketed  factor  may 
be  approximately  unity.  Then  apart  from  the  basic  absorption  matrix  element, 
(and  Planck's  constant),  one  has 


10  v,  a(u>-w  ) 
in  o 


(1.72) 


It  will  now  be  shown  that  the  r.h.s.  of  (1.72)  is  equal  to  the  time  during  which 
the  collision  produces  a  momentary  vibrational  frequency  perturbation,  which 
lies  within  a  unit  differential  range  about  the  spectral  shift,  w-wo.  One  notes, 
from  (1.16)  (and  Bohr's  frequency  condition)  that 


«  [V11(R)  -  Voo(Il>] 


(1.73) 


which,  with  use  of  (1.28a)  and  (1.28b),  becomes 


,  ..  -aR 

0)p  -  (u^-u^)  e 


(1.74) 


(1.75) 


A  "frequency-coincidence" 


a)  (R)  “  u>— u)  , 

p  o’ 


is  achieved  at  a  value  of  R  given  by 


d  I  i  VMi 

a)  a  °®  ai-w 

o 


(1.76) 
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The  differential  range,  dR^  in  which  cu  (R)  will  lie  with  a  differential  frequency- 
range,  dto  is  clearly 

dR  ■  da)/a(a)-ci)  )  ;  (1  •  77) 

U!  O 

correspondingly,  the  time,  dt,  during  which  ^(R)  will  be  contained  within  a  unit 
frequency-range  about  is 


dt 

dot 


in 


dR 

u 

du 


v  o(u>-u>  ) 
in  o 


(1.78) 


q.e.d.  The  coincidence  of  (1.72),  and  hence  (1.71),  with  the  prediction  of  the 
statistical  theory  is  thus  established. 


Further  detailed  discussion  of  the  statistical  theory  will  be  abandoned 
for  the  time  being.  The  motivation  for  this  decision  is  that,  as  will  now  be 
shown,  it  appears  that,  under  most  circumstances  the  frequency  perturbations 
required  for  the  far-wing  absorption  of  Interest  in  these  studies  (600-1000  cm  *) 
are  simply  not  attainable. 

The  consideration  begins  by  noting  that  [as  given,  e.g.,  by  (1.66)], 

2 

k  be  real,  i.e.,  <  >0.  Applying  this  condition  to  the  r.h.s.  of  (1.66)  one 

readily  obtains  the  inequality 


k*  -n(u)-u)  )+E. 
f  f  o  i 


or,  from  (1.40) 


fi(a)-a)  )+E. 
o  l 


l+3a2xV/4y  u 

V  o 

1+a^X^h/Ay  a > 
v  o 


(1.79) 


From  (1.79),  it  is  clear  that  for  the  determination  of  the  domain  of  applicability 
of  the  statistical  theory,  the  value  of  the  parameter 


.  _  aV* 

5  *  ’ 

V  o 

is  of  decisive  importance. 


(1.80) 
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In  order  to  estimate  £,  let  us  consider  HC1  as  a  representative  diatomic 
molecule.  For  this  case,  if  it  is  considered  that  atoms  A  and  B  (referred  to  in 
e.g.  (1,4)  and  following  text)  are  respectively  Cl  and  11,  one  may  take  X  =  1. 

Moreover,  u  may  be  identical  with  the  mass  of  the  hydrogen  atom;  thus 

-24  13 

u  =  1.6x10  g.  For  ui  ,  Herzberg's  tables  yield 
v  o 

3  -1 

a)  =  2itc  3x10  sec 
o 


Finally,  for  a,  let  us  use,  as  a  representative  value,  the  mean  of  those 
quoted  in  Chapter  3,  page  71  of  this  report  (a  is  denoted  there  by  the  symbol  X/2) 

’  /  in8  -1 

a  =  4x10  cm 

One  then  obtains 

£  **  1/12  .  (1-81) 

Introducing  this  value  into  (1.79)  yields 

h(u)-w  )  ;  (1.82) 

o 

taking  E  ~  kt  —  300  cm-1,  it  is  clear  that  one  obtains  for  the  frequency  range 
i  _  ^ 

of  applicability  of  the  statistical  theory,  a  value  ~  30  cm  . 

Without  proceeding  further,  it  should  nevertheless  be  pointed  out  that  the 
contributions  of  superthermal  incident  energies  should  not  be  dismissed  out-of¬ 
hand.  Although  adversely  weighted  by  the  Boltzmann  factor,  exp (-E^/k't)  ,  they  may 
still  make  a  respectable  showing  in  comparison  with  the  other  (thermal)  con¬ 
tributions,  which  fall  outside  the  range  of  the  statistical  theory;  as  shown  in 
the  following  paragraphs,  these  are  also  of  exponential  smallness  in  the  fve- 
quency  shift. 

After  these  remarks  let  us  explore  the  case  in  which  [referring  to 
Eq.  (1.44)]  there  are  no  saddle  points  on  the  real  axis  (from  x=l  to  x=«>). 
Physically  this  case  corresponds  to  the  situation  described  at  the  end  of  section 
I,  in  which  the  collision-induced  vibrational-frequency  perturbation  is  too  small 
to  provide  the  momentary  coincidence  with  the  spectral  shift,  u)-u)q,  which  is 
required  by  the  statistical  theory. 
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For  the  case  of  saddle  points  off  the  real  axis,  it  can  be  shown14  that 
the  integration  contour  can  be  deformed  so  that  the  principal  contribution  to 
the  integral  occurs  in  the  vicinity  of  the  saddle  point  in  the  upper  half  of  the 
complex  plain;  from  Eq.  (1.55)  this  saddle  point  is  given  by 

zs  "  fU+1/2)  +  i[c/Y-(c2/4)(l+l/Y)2]1/2  .  (le83) 


The  calculation  of  the  saddle-point  integral  proceeds  quite  analogously 
to  the  case  of  real  saddle  points.  The  final  result  is  [compare  with  (1.62)] 


Z 


j  v+1 

1  6  (u+v)tt/2 


sin  J(Tr/ 2)  (p+v)J 
where  [cf.  (1.51)  and  (1.83)] 


exp 


2r-, 


1 1/2 


Mu+v) 


(1.84) 


-  -1^1  f  -  -(kf+ki)f/a 


(1.85a) 


tan 


tao'1  [t<T/y-1)]+  * “"'‘[.(Wrt-i]  • 


(1.85b) 


and 


6  "  [C/Y  -  (c2/4)(1+1/y)2]  /  • 


(1.85c) 


From  (1.50),  (1.48a, b),  (1.40),  and  (1.59b)  one  obtains  [by  algebraic  manipulations 
similar  to  those  which  led  to  (1.65)] 


V  1  Vf 
*  2  ..  2 


tt3/4 


<Vf-Vi)  sinh^(ir/a)(kf+k1) 


k2_k2^ 

f  i 


<o 


exp[-2(kf+  k±)f/a]  , 


(1.86) 


where 
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kf~kl  2 
1_Vi/Vf  "  f 
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1-vi/vf 


1+E  ) 
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2M(ai-to  )  f 

- — — —  1/C  "  E  /li (to— (jj  )  +  ■=• 

n  1  O  z 


(1.87) 


with  5  being  given  by  (1.80),  i.e.. 


5  “  2y  u  5 

v  o 

as  pointed  out  in  the  text  subsequent  to  Eq.  (1.80),  its  numerical  value  may 
generally  be  expected  to  be  small  compared  to  unity. 

Upon  substituting  (1.86)  into  (1.38),  one  has 
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r  £uo  1 
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TT 

f2M(oi-uio)l 

i — * 
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3 

w 

)C 

a 

[  hn  | 

exp  (-  [8M(ui-wo)  ^  g)  , 


(1.88) 


where  the  wiggle  above  the  equality  sign  denotes  an  approximation  therein  a  factor 
equal  to  the  curly  bracket  of  (1.71)  has  (in  accordance  with  the  text  subsequent 
to  that  equation)  been  replaced  by  unity, 


n  =  1/5  -  E1/ft(u-uQ)  +  1/2, 


[with  5  being  given  by  (1.80)  and  (1.81)],  and 


r, 


E  /K  (a)-(i)  ) 

i  o 


1/2 


1+E  /fi  (ui-m  ) 
i  o  } 


1/2 


(1.89) 


(1.90) 


It  is  now  of  interest  to  "normalize"  the  r.h.s.  of  (1.88)  with  respect  to 
the  standard  Lorentz-type  absorption  profile.  To  this  end,  let  us  briefly  sketch 
a  derivation  of  the  phenomenological  impact-theory  formula  for 
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One  considers  the  diatomic  molecular  oscillator  (initially  in  its  ground 
state)  in  the  presence  of  the  electromagnetic  field  of  the  light  wave  and  subject 
to  certain  collisional  events  which  have  the  following  two  properties : 

1)  They  are  random  in  time;  that  is  to  say,  a  collision  occurs  in 
a  differential  interval,  dt,  at  a  time  t  after  the  preceding 
collision  with  the  probability 


dW 


coll 


e-t/T dt/T 


(1.91) 


2)  Their  effect  is  to  randomize  the  phase  relationships  between 
the  ground  and  excited  states  of  the  vibrator.  As  is  known, 
such  phase  relationships  occur  as  a  result  of  the  perturbative 
action  of  the  electromagnetic  field;  in  particular,  ac  a  time, 
t,  after  the  electromagnetic  field  is  "turned  on"  (i.e.,  after 
a  previous  randomizing  collision),  the  coherent  part  of  the 
vibrational  wave  function  is^ 


*o  + 


^W10 

-iwt  -iu  t)  M10  f  +iwt  -iw-tl 

15(u)-00  ) 

0 

,e  "e  °  j+-K(o,-h,o)  (e  ~e  °  j 

(1.92) 


The  significance  of  property  (2)  is  that,  immediately  after  a  collision 
the  coefficient  of  <j>  in  (1.92)  gets  multiplied  by  an  arbitrary  phase  factor.  It 
then  follows  that  the  state  of  the  system  is  specified  entirely  by  the  occupation 
probability  of  the  individual  eigenstates,  in  particular,  if  the  collision  occurs 
at  time  t  (subsequent  to  the  previous  one),  the  occupation  probability  of  the 
first  vibrational  level  is  augmented  by  an  amount  equal  to  the  absolute  square 
of  the  coefficient  of  in  (1.92).  If  this  quantity  be  averaged  over  the  time, 
t,  between  collisions,  with  a  weight  factor  given  by  the  r.h.s.  of-  (1.91),  one 
obtains  the  average  probability  per  collision  that,  as  a  result  of  the  simul¬ 
taneous  action  of  collisions  and  the  electromagnetic  field,  a  (0  1)  vibrational 

transition  occurs.  Carrying  out  the  indicated  operation,  one  obtains  for  this 
probability  a  sum  of  three  contributions,  two  of  which  are  proportional  to  the 
absolute  squares  of  the  first  and  second  term  of  the  square  bracket  of  (1.92), 
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respectively,  the  third  being  the  interference  fern.  In  what  follows,  let^us 
restrict  the  discussion  to  the  "resonance"  region  defined  by  the  condition 


In  this  case,  the  first  contribution  dominates,  thereby  permitting  us  to  consider 
the  transition  as  being  associated  with  electromagnetic  absorption.  In  conformance 
with  this  feature,  let  u«i  denote  it  as  (the  superscript  "(im)"  referring  to 

the  impact-theory  basis  of  its  derivation). 

In  carrying  out  the  average  over  t,  let  us  introduce  the  further  condition 


L-U)  I  >>  1/t  , 

'o' 

(which  is  appropriate  for  the  discussion  of  the  wings  of  a  Lorentzian  line,  and 
which,  moreover,  is  quite  well  obeyed  in  atmospheric  environments  for  the  spectral 
shifts  of  relevance  to  the  present  paper).  One  then  obtains 


(im) 

10 
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|1I(u>-u)o)J 


2 


(1.94) 


It  is  now  necessary  to  recognize  that  the  quantity  of  ultimate  interest  is 
not  the  absorption  probability  per  collision,  but  rather  the  absorptive  proba¬ 
bility  per  unit  time  (which  will  here  be  denoted  as  W1Q) .  Generally,  the  relation 
between  the  two  quantities  is  of  the  form 


where  v  is  a  suitably  defined  "effective"  collision  rate.  In  the  case  of  the 
phenomenological  impact  treatment  sketched  above,  vgff  is  simply  the  reciprocal 
of  t,  so  that 
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From  (1.96),  It  is  now  apparent  how  to  rewrite  the  basic  result,  (1.88), 
of  the  present  calculation  in  a  form  suitable  for  comparison  with  the  impact 
theory;  namely, 


(1.97) 


(1.98) 


where  is  the  effective  number  of  collisions  per  second. 

In  a  one-dimensional  model  of  relative  translational  motion  —  which,  in 
effect,  has  been  imposed  on  us  by  the  collinear  collision  restriction  —  v 

ef  f 

would  be  something  like  the  relative  velocity  of  molecule  and  perturbing  atom 
times  a  suitably  defined  linear  density  of  perturbing  atoms.  In  order  to  go 
behond  a  phenomenological  prescription  of  this  type,  it  is  necessary  to  generalize 
the  treatment  to  include  collisions  other  than  collinear.  A  crude  preliminary 
treatment,  which  should  (in  the  opinion  of  the  present  author)  nevertheless  yield 
results  of  semiquantitative  significance,  will  now  be  given. 

Specifically,  it  will  still  be  assumed  that  the  trajectory  of  the  colliding 
atom,  relative  to  the  center-of-gravity  of  the  diatomic  molecule  is  parallel  to 
the  diatomic  axis  —  the  direction  of  said  axis  being  taken  as  fixed,  in  accordance 
with  the  neglect  of  rotational  motion  (which  was  introduced  at  the  beginning  of 
the  treatment).  However,  the  distance  between  the  trajectory  and  the  center-of- 
gravity  (that  is,  the  classical  impact  parameter)  is  now  permitted  to  take  on  a 
nonvanishing  value. 

A  simple  procedure  for  the  treatment  of  these  "off-center"  collisions  is 
afforded  by  an  approximation  —  here  designated  as  the  Takayanagi  approximation, 
in  view  of  its  prior  use  by  that  author  in  a  paper  on  vibrational  deactivation^ 
which  is  nevertheless  (in  the  opinion  of  the  present  author)  quite  suitable 
for  obtaining  results  of  semiquantitative  accuracy.  In  describing  the  approxi¬ 
mation,  let  us  first  note  that,  for  impact  parameters  t  atomic  radii,  (~3X),  the 
incoming  particle  finds  itself  in  a  spherical  potential  (centered,  of  course, 
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around  atom  B  of  the  molecule).  It  is  to  be  noted  that  the  scale  of  variation  of 
the  potential  is  ~l/a  25&.  This  dimension  is  substantially  smaller  than  the 
distance  of  closest  approach,  (~3&).  In  other  words,  the  radius  of  curvature  of 
the  equipotential  surfaces  is  large  compared  to  the  distance  (~l/a)  in  which  the 
colliding  particle  experiences  its  principal  interaction  with  the  molecule.  In 
the  Takayanagi  approximation  the  curvature  of  these  surfaces  is  neglected.  Speci¬ 
fically,  the  potential  is  taken  as  one-dimensional  with  the  equipotential  planes 
oriented  perpendicular  to  that  radius  which  intersects  the  trajectory  at  its 
distance  of  closest  (classical)  approach.  The  potential  on  each  plane  is  taken 
to  be  that  of  the  (actually)  equipotential  sphere  which  is  tangent  to  the  plane  in 
question.  With  this  approximation,  one  writes  the  wave  function  as  a  product  of 
two  factors,  describing  motion  parallel  and  perpendicular  to  the  (fictitiously) 
equipotential  planes,  respectively.  The  parallel  motion  is  described  by  a  plane 
wave  whose  wave  vector  lies  in  the  "plane  of  incidence"  (i.e.,  the  one  containing 
the  classical  trajectory  and  the  center  of  atom  B)  and  has  a  magnitude 


kM  *  k8in9  "  kP/*c  »  (1.99) 

(where  0  is  the  angle  between  the  trajectory  direction  and  the  radius  which 
contacts  it  at  the  distance  of  closest  approach,  rc,  and  p  is  the  distance  from 
the  trajectory  to  the  center  of  atom  B) .  (For  the  particular  collisions  under 

consideration,  p  is  also  the  distance  from  the  trajectory  to  the  molecular  center 
of  gravity) . 


For  the  perpendicular  motion,  one  is  then  left  with  a  SchrtJdinger  equation 

which  is  identical  to  (1.31)  except  that  k  is  to  be  replaced  by  its  "perpendicular" 
component'1' 


k  cos  0 


l-p2/r2l 


1/2 


ci 


(1.100) 


Equivalently,  the  incident  energy  Ei  undergoes  the  replacement19 


Ei  + 


cos2  0 


(1.101) 
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From  the  above  results  [In  particular,  from  (1.01)  and  pnceding  text],  it 
follows  that  the  transition  probability  at  an  Impact  parameter  p  obeys  the 
equation 


P10(P’V 


10 


0,E  (l-p2/r2) 
1  c 


(1.102) 


With  the  aid  of  (1.102),  the  approximate  first-principles  replacement  of 
the  phenomenological  equation  (1.98)  reads 
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N  TTV, 


2a 


2irpdp 


(1.103) 


where  Np  is  the  density  of  perturbing  atoms,  =  (2El/M)1/2  is  the  initial 
relative  velocity  and  where  n(o)  and  g(p)  are  gotten  by  using  the  replacement 
(1.101)  in  (1.89)  ai.d  (1.90)  respectively. 

A  quantitatively  accurate  evaluation  of  (1.103)  for  all  cases  of  interest 
remains  to  be  carried  out  in  the  future.  Preliminary  studies  indicate  that,  if 
(1.103)  be  expressed  in  the  form 
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8  »P2 
opt 


(1.104) 


[where  n  and  g  are,  as  before,  the  values  of  n(p)  and  g(p)  for  p  -  0,  and  where 

Popt’  the  °Ptlcal"  collision  radius,  is  of  course  to  be  determined  by  integrating 
(1.103)],  then 


Popt  *  rco  ~  '  (1.105) 

That  is,  the  upper  limit  of  the  optical  collision  radius  is  of  the  order  of  a 

gas-kinetic  value.  It  will  now  be  shown  that,  even  with  this  upper  limit,  the 

far  wing  broadening,  given  by  (1.104)  is  numerically  far  below  the  impact-theory 
result. 
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In  carrying  out  the  numerical  comparison,  let  us  approximate  q  by  its  first 
term,  which  is,  cf.(1.81),  large  compared  to  unity  (the  other  terms  being  numer¬ 
ically  *1);  in  particular,  let  us  fix  n  at  the  value  suggested  by  (1.81),  i.e., 

*  "  10  •  (1.106) 


Then.,  upon  using  the  estimate 

„  .  in8  -1 
a  ■  4  x  10  cm  , 

and  the  relationship 

(2E  )  r  vl/2 

vi  "  M  *  M  -  2.2  Xl05*t1/2, 

(where  =  - -  and  where  T  has  been  set  equal  to  300°K) , 

proton 

Np  -  2.7  x  1019^  cm"3  , 

(where  7l  is  the  perturbing  gas  density  in  atmospheres),  one  has 


(1.107) 


(1.108) 


as  well  as 


(1.109) 
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(1.110) 


where  Av  is  the  frequency  shift,  u>— u>o ,  expressed  in  wave-numbers,  and  where 

is  given  in  angstroms;  in  what  follows  p  _  is  set  equal  to  three, 
opt  opt 

Turning  to  the  evaluation  of  gr  as  given  by  (1.90)  and  (1.85b),  one  easily 
verifies  that,  with  the  ^-lantity  c  ■  V^/V^  “  K  set  at  a  fixed  value,  (1/10),  g 
is  only  a  function  of  the  dimensionless  parameter,  ft(w-u>o)/E^.  For  a  repre¬ 
sentative  value  of  this  parameter,  say  3,  (corresponding  to  Av  900  cm  ^) ,  a 
numerical  calculation  gives  g  ■  1.1.  One  then  has,  even  for  the  extreme  case  of 
M  ■  1  (mass  of  relative  motion  ■  hydrogen  mass) 

e-4»l1/2(M;/103)1/28  ^  e-4  .  _018  _ 
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The  corresponding  value  of  1/t  is  found  to  be 


T  ■  (*9)  (1.8)  x  10  (9)(.018)  ■  2.8  x  108sec  3  . 

Expressed  in  wave-number  units,  this  is  2.8  x  108/2irc  *  1.4  x  10~3cnf3,  a 
substantially  below  the  corresponding  quantity  used  in  Penner's  article,3 
4  x  10  cm  ^/atmosphere. 

Another  ..representative  calculation:  for  K(aj-w  )/E  =  2(Av  •  600  cm~ 
one  has  g  -  .88;  exp^^/^Av/lO-^)1/^]  -  exp[-4(?6)li/2(. 88)]  -  e‘2*72 

The  corresponding  value  of  1/t  is 

~  "  (*6)  (1.8)  x  10  (9) (.066)  ■  8  x  108sec  3  ■  4  x  10-3cm""3 

—  a  factor  of  ten  less  than  the  Penner  value.23 


value 
namely , 


-  .066. 
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SUMMARY  AND  FUTURE  PROSPECTS 


In  this  paper  a  model-type  treatment  of  the  far  wing  broadening  of  the 

vibrational  lines  of  a  diatomic  molecule  by  a  monatomic  perturber  has  been  given. 

The  model  is  one  in  which  interatomic  forces  are  assumed  to  be  of  the  exponential 
— ar 

type,  e  ,  and  in  which  rotational  motion  is  ignored.  Specific  results  are 
derived  for  collinear  collisions;  these  are  generalized  somewhat  by  the  use  of  an 
approximation  due  to  Takayanagi. 

It  is  found  that  the  so-called  statistical  theory  is  valid  only  out  to 
distances  of  the  order  of  30  cm  1  from  the  line  center.  Beyond  this  point  the  line 
profile  falls  exponentially  with  increasing  frequency  shift.  Numerical  estimates 
at  distances  —'600  cm  and  900  cm  ^  from  the  line  center  indicate  that  the  actual 
absorption  falls  substantially  below  the  impact-theory  estimates  of  Penner.^ 

With  regard  to  future  work  in  this  area,  the  following  items  may  be  listed: 

1)  First  of  all,  a  systematic  numerical  analysis  of  the  final  formulae 
of  the  present  paper  is  required,  in  order  to  present  a  proper 
picture  —  in  particular  the  overall  frequency  variation  of  the 
line  profile,  and  a  comparison  with  the  impact-theory  predictions 
of  the  model  (exponential  interatomic  forces) . 

2)  The  role  of  attractive  forces  should  be  considered;  it  may  be 

possible  to  do  this  rather  straightforwardly,  using  Morse-type 
22 

potentials.  Quantitatively,  attractive  forces,  by  speeding  up 
the  particles,  increase  the  effective  incident  energy;  this  increase 
is  preliminarily  expected  to  lead  to  an  enhancement  in  the  ab¬ 
sorption  probability. 

3)  By  an  extension  of  the  Takayanagi  approximation,  it  appears  pos¬ 
sible  to  generalize  the  treatment  further,  to  include  collisions 
where  trajectories  are  inclined  at  an  arbitrary  angle  to  the 
molecular  axis. 

4)  It  would,  of  course,  be  eminently  desirable  to  take  account  of 
rotational  motion.  This,  however,  promises  to  be  a  tricky  business. 
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3.  These  magnitudes  are  suggested  in  the  article  of  S.  Penner,  Laser  Summer 
Study,  (1965),  The  Institute  for  Radiation  Physics  and  Aerodynamics, 

Univ.  of  Calif,  at  San  Diego,  in  connection  with  the  absorption  of 

1  micron  due  to  H^O  lines  situated  in  the  spectral  region. 

4.  It  may  at  this  point  be  mentioned  that  realistic  value  for  a  lie  in  the 

8  —1 

range  4  -  5  x  10  cm  .  Such  values  begin  to  approach  the  requirements 
of  condition  (a) . 

5.  cf.  e.g.  G.  Herzberg,  Spectra  of  Diatomic  Molecules  (Van  Nostrand,  1950), 
pp.  125-127,  especially  the  quoted  values  of  the  frequency  separations 

of  the  maxima  of  the  two  rotational  branches  (P  and  R)  for  HC1  (124  cm  ^) 
and  CO  (55  cm  ^)  at  room  temperature.  These  values,  of  course, 
characterize  the  rotational  states  of  the  molecule  before  collision; 
after  collision,  the  "effective"  rotational  temperature  will  in  general 
be  augmented  by  the  conversion  of  the  energy  K(u>-u>o)  into  rotational 
plus  translational  energy.  However,  since  room  temperature  corresponds 
to  ^300  cm  \  the  values  quoted  here  will  not  be  augmented  by  more  than 
a  factor  of 

((u>-wo)  +  kT 
kT 

6.  Alternatively  stated,  the  conversion  of  an  appreciable  fraction  of  the 
energy  (w-wo)  into  energy  of  rotation  would  require  AJ  >>  1. 
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7. 


::r;“ hara  that  in  une  “uh  thi°  «-  «»»i 

8  '  °  r°tatlonal  factions  are  easily  shewn  to  be  comparable 

Irialn  tI“1“1“l  f“8’  b°th  *“  — —  -  energy 

cf.  N.F.  Mott  and  8.W.  Massey,  Theory  of  Atomic  Collisions,  3rd  Ed 
(Oxford  1965)  pp.  686-690. 

It  may  parenthetically  he  remarhed  that  the  approximation  represented 

y  qs.  1.8  1.10)  -  the  so-called  distorted-wave  approximation  - 

should  give  good  results  since  (a)  in  practice,  by  virtue  of  the  repulsive 

Character  of  the  interaction,  V(R,X)  will  generally  he  of  the  order  of 

the  incident  energy  E,  u  kT  «  Bw,  (larger  values  of  the  repulsive 

potential  being,  in  essence,  dynamically  inaccessible).  As  will  later 

be  seen  the  intervibrational  matrix  elements  of  V(R,»  .re  even  smaller. 

s  t  e  neglect  of  these  matrix  elements  which  constitute  the  distorted 
wave  approximation . 

Here,  it  is  necessary  to  introduce  the  full  time  dependence  of  the 
unperturbed  wave  function. 

and'th11  TaUy  ^  tha  ““  th8t  the  °£  «■  range  parameter,  a, 

then  ‘  ^“'ion-amplitude  scale-factor,  «./„„«,/«  l8  .mall.  One  may 

expand  the  exponential  factor,  l«»»,  i„  the  integrals,  obtaining 
Vi  ■  Vo  “»<•»».>  [hoVb/A^J,  vf  .  Vo  exp(„AYc)  [iWoVb/^wJ. 

th'  8rPtlrn8ly’  the  r-h'8-  °£  (1-76>  is  identical  with  the  R  for  which 
the  above  discussed  modification  of  the  Franok-Condon  Principle  Tcf 

Eqs.  (1. 67a,b)  and  (1.68)1  is  satisf^H  u 

n  18  satisfied  (as  may  be  easily  verified  by 
the  reader) .  J 

o.  Herrberg,  Spectra  of  Diatomic  Molecules  (Van  Nostrand),  2nd  Ed.  (1,50) 

p.  534.  '  * 

The  details  will  have  to  be  relegated  to  a  future  report. 

The  initial  time  dependence  e~^^it/R 

,H11  .  P  ’  S  *  occurs  as  a  common  factor,  and 

•*ill  hence  be  discarded. 
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16.  It  may  be  remarked  in  passing  that,  if  this  condition  is  not  fulfilled, 
the  randomization  recipe  of  (2)  is  not  correct.  Namely,  as  shown  by 
Van  Vleck  and  Weisskopf  [Rev.  Mod.  Phys.  1J,  227  (1945)],  the  total 
vibrational  state  immediately  after  collision  is  to  be  considered  as  an 
incoherent  aggregate,  not  of  the  eigenstate  of  the  field-free  vibrator, 
but  rather  of  the  vibrational  eigenstates  in  the  presence  of  a  static 
electric  field,  equal  to  the  momentary  value  of  the  oscillating  electro¬ 
magnetic  field  at  the  impact-time,  t.  A  detailed  analysis  shows  that 
this  recipe  becomes  equivalent  to  the  one  used  above  when  (1.93)  is 
fulfilled. 

17.  K.  Takayanagi,  Prog.  Theoret.  Phys.  497  (1952). 

18.  There  is,  of  course,  a  further  correction  because  of  the  fact  that  the 
direction  of  molecular  vibration  no  longer  coincides  with  that  of  the 
line-of-centers  of  the  colliding  atoms  (B  and  C).  This  involves  alter¬ 
ations  in  the  formulae  of  Ref.  11,  for  the  constants  and  V^.  These 
alterations  turn  out  to  be  of  secondary  importance  and  will  be  ignored 
in  this  report. 

19.  It  may  here  be  remarked  that  Takayanagi 's  formulation  is  somewhat 

different  than  the  foregoing.  He  separates  the  Schrddinger  equation  of 

relative  motion  in  spherical  coordinates,  obtaining  a  radial  equation 

which  differs  from  (1.31)  in  the  presence  of  the  usual  centrifugal  term, 

2 

fc()l+l)/r  .  Takayanagi' s  approximation  is  the  replacement  of  this  term  by 
2  2  2 

£(£+l)/rc  pt  1  /rc,  which,  with  the  introduction  of  the  impact  parameter, 
p  = 'K£,/mv,  is  easily  seen  to  be  equivalent  to  (1.100). 

20.  It  is  incidentally  to  be  noted  that  r£  is  itself  a  function  of  p,  to  be 
determined  by  solution  of  the  equation  [cf.  (1.31),  (1.30c),  and  (1.1)] 

0  =  k2(l-p2/r2)  -  (2M/fi2)vo  e_0trc  • 

However,  for  the  case  in  which  p  <<  r  (where  r  is  the  value  of  r 

co  co  c 

for  p  =  0) ,  one  may  clearly  approximate  r£  by  rco<  As  will  shortly  be 
seen  below,  only  this  case  is  treated  explicitly  in  the  present  paper. 


It  should  be  kept  in  mind  that  the  choice^-  1  yields  an  extreme  upper 
limit  for  1/t.  E.g.,  for  HjO,  it  would  seem  that  even  with  the  rotational 

motion  taken  into  account,  a  minimum  value  of#l“  2  would  seem  to  be  more 
appropriate.  For  this  choice,  the  value  for  1/t  is  reduced  by  a  factor 
e-2.72(.41)  x  33>  The  resuit9nt  vidth  is  30  times  smaller  than  the 

Penner  value. 

Cf.  Devonshire,  Proc.  Roy.  Soc.  A158 ,  269  (1937). 
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Chapter  2. 


ABSORPTION  STUDIES  OF  THE  (0^  COMPLEX 


( 


This  chapter  describes  studies  of  various  physical  processes  involved  in  tne 
heating  of  the  atmosphere  in  the  vicinity  of  a  high-powered  laser  pulse.  The 
starting  point  of  these  studies  is  to  be  found  in  Appendices  E  and  G  of  the 
~—er  — rm^er-  Stiidy  Report  (K.A.  Brueckner,  Institute  for  Radiation  Physics  and 
Aerodynamics,  UCSD,  La  Jolla,  California,  Aug.  9-20,  1965).  The  analyses  reported 
in  the  latter  assume  as  an  engineering  requirement  that  the  effective1  absorption 
coefficient  be  less  than  or  of  the  order  of  10"llcm-1  (at  pressures  of  the  order 
of  one  atmosphere).  Specifically,  Penner  and  Olfe  (Appendix  E,  loc.  cit.)  arrived 
at  tha  conclusion  that  in  the  vicinity  of  the  wavelength  10,468  X  the  absorption 
coefficient  k(v)  is  ~1.5  x  10  11cm  1.  With  this  result  the  engineering  problem 
would  appear  to  have  been  solved,  were  it  not  for  the  fact  that,  as  pointed  out  by 
Kolb  and  Ory  (cf.  Appendix  G,  loc.  cit.),  atmospheric  absorption  in  the  spectral 
neighborhood  of  10,468  X  is  actually  dominated  by  a  continuous  absorption  band, 
due  to  the  existence  (momentary  or  otherwise)  of  (02)2  complexes.  From  the  data 
of  Dianov-Klokov,  Kolb  and  Ory  conclude  that  at  10,468  X  the  absorption  coef- 
ricient  is  1.5  x  10  cm  1  (down  by  a  factor  of  ten  from  the  peak  at  10,600  X). 

It  is  clear  that,  if  the  energy  absorbed  from  the  laser  beam  is  immediately 
converted  into  heat,  one  is  reduced  to  the  investigation  of  (a)  the  possibilities 
of  an  atmospheric  window  in  the  further  wings  of  the  (02)2  band  (e.g.,  in  the. 
neighborhood  of  say  10,300  X  or  10,800  X),  or  (b)  looking%or  a  window  in  an 
entirely  different  spectral  region  (such  as  the  10  micron  window  mentioned  by 
Kolb  and  Ory).  Deferring  these  approaches  for  the  time  being,  let  us  consider 
in  some  detail  the  physical  processes  responsible  for  the  conversion  (wholly  or 
in  part)  of  the  absorbed  laser  energy  into  heat. 

To  begin  with  one  may  note  from  the  results  of  Badger  et  al4  (as  well  as 
from  the  earlier  references  quoted  therein)  that  the  10,600  X  band  is  associated 
with  the  transitions 

(\«».  3yo)j  *  (3yo),  ^(ij) 

pyi).  x*.(0>)  . 
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in  which  the  resulting  complex  may  be  regarded  as  consisting  of  a  pair  of  loosely 

bound  (or  momentarily  colliding)  diatomic  oxygen  molecules.  One  of  these  molecules 

ends  up  in  the  first  excited  electronic  state  ,  the  other  in  the  electronic 
3  8 

ground  state  Z  ;  in  addition,  one  of  the  molecules  acquires  a  single  quantum  of 

O 

vibrational  excitation.  As  will  be  seen  shortly,  there  is  good  evidence  that  under 
atmospheric  conditions  the  lifetime  of  the  electronic  state  is  considerably  in 
excess  of  one  millisecond;  on  the  other  hand,  there  is  equally  good  evidence  that 
(under  typical  mid-latitude  atmospheric  conditions)  vibrational  de-excitation  of 
O2  takes  place  in  a  time  substantially  less  than  one  millisecond. 

Before  getting  into  the  detailed  estimation  of  relaxation  times,  let  us 

briefly  note  some  interesting  (and  perhaps  eventually  pertinent)  facts  concerning 

the  structure  of  the  (02)2  dimers.  According  to  Ref.  4  [cf.  their  Eq.  (6)]  the 

radiative  decay  rate  of  ^A  states  is  given  by  the  formula 

S 


-  (sec-1) 


■  2.6  x  10 


-4 


1  +  3.8  PA  +  J.O  P„„  +  .7  P 


0 


CO, 


N, 


where  the  P's  are  partial  pressures  in  atmospheres.  Of  interest  here  is  the  term 

proportional  to  P«  ;  it  gives  the  strength  of  the  radiative  transition  associated 
u2 

s;ith  the  (0^2  complex  (either  a  stabilized  dimer,  or  a  momentary  0^  complex 
formed  by  two  colliding  0„  molecules).  Now,  the  likelihood  of  an  09  molecule 
being  in  the  atomic  vicinity  RQ  _Q  -2  x  10“  cm  is  ~  NQ  R^  _q  .  Introducing 
the  symbol  1/Tq  to  denote  an  effective  radiation  rate  for^the^  ^ 


0^  complex, 


we  then  have 


(2.6  x  10  4) (3.8)  P 


0, 


1  hr 


Ai  n  R3 
3  N02  R02-02. 
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With  Nq  3  x  10  Pn  ,  and  with  the  above  estimate  for  Rn  „  one  then  has 


■0 


°2-°2 


1/t  ~  1  sec  1, 


which  is  extremely  small  compared  to  the  rate  of  allowed  electric  dipole  tran¬ 
sitions.  This  result  suggests  an  extremely  loose  coupling  of  the  two  02  molecules, 
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such  that  the  spins,  parities,  and/or  axial  angular  momenta  of  the  individual 
diatomic  constituents  remain  good  quantum  numbers. 

This  conclusion  is  supported  by  the  very  low  quoted  estimates  of  the  dis¬ 
sociation  energies  of  ( dimer  states.  According  to  Bader  and  Ogryzlo  these 
range  from  200-600  cal/mole,  which  is  equivalent  to  .009-. 025  ev.  Such  small 
dissociation  energies  are  indicative  of  lo;  g  range,  extremely  weak,  attractive 
interactions  (probably  of  van  der  Waals  origin)  in  which  the  individual  02  mole¬ 
cules  are  kept  well  apart  by  the  usual  shorter  range  but  relatively  strong 
repulsive  interaction. 

It  may  be  remarked  here  that  Bader  and  Ogryzlo  emphasize  that  their  results 
lead  to  the  conclusion  that  the  radiating  complexes  are  true  dimers.  However,  in 
our  opinion,  the  low  dimer-dissociation  energies  quoted  by  them  strongly  suggest 
that  such  dimer  states  exist  for  at  most  a  few  collision  times  '*’10  sec.  For 
the  overwhelming  fraction  of  time,  electronic  and  vibrational  excitations  reside 
in  individual  02  molecules. 

With  these  preliminary  observations  out  of  the  way,  let  us  proceed  to  a 
discussion  of  the  lifetime  of  the  electronic  state  1Ag  under  atmospheric  conditions. 
Here,  very  fortunately,  relatively  strong  experimental  evidence  for  a  lifetime  far 
in  excess  of  one  millisecond  is  provided  by  the  previously  quoted  work  of  Bader 
and  Ogryzlo.  These  authors,  using  calorimetric  techniques,  find  (cf.  curve  (b)  of 
Figure  1)  that  in  a  fast-flow  system  in  which  gas  flows  from  a  discharge  at  a 
drift  speed  vd  ~  200  cm/sec,7  the  concentration  of  excited  02  molecules  (which 
they  identify  as  XA  molecules)  does  not  exhibit  perceptible  diminution  in  a  flow- 
distance  of  45  cm;  moreover,  said  persistence  of  excitation  is  unaffected  by  the 
presence  of  small  amounts  (~1  percent)  of  H20,  N02,  N20  (as  well  as  some  non- 
atmospheric  constituents).  One  may  immediately  infer  that  under  these  experimental 

conditions 


Postulating  a  hypothetical  volume-destruction  process  involving  the  binary  col 
lision  of  excited  02  molecules  with  other  atmospheric  constituents  (02  or  N2) 
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one  may  then  infer  a  lifetime 


>>  .23(4.4)/760  =  1.3  millisecond. 

Ag 

Two  possible  loopholes  in  the  above  chain  of  reasoning  may  be  pointed  out 
here.  One  is  the  tacit  assumption  that  0*  -  N2  collisions  will  not  be  more  ef¬ 
fective  than  0*  -  02  collisions  for  electronic  de-excitation.  In  our  opinion  this 
assumption  appears  [in  the  light  of  preliminary  study  of  the  theory  of  such 
quenching  reactions,  as  presented,  e.g.,  in  the  paper  by  K.  Laidler,  J.  them.  Phys 
34  (1942 (]  to  be  eminently  reasonable.  Similarly,  it  appears  quite  reasonable 
to  assume  that  minor  atmospheric  constituents,  such  as  carbon  dioxide,  will  not  be 
more  effective  than  H20  in  excitation-quenching. 

The  other  loophole  is  the  assumption  that  the  volume-destruction  process  is 
two-body.  If,  on  the  contrary,  it  were  a  three-body  process,  one  would  have  to 
multiply  the  above  estimate  of  1.3  milliseconds  by  another  factor  of  4.4/760, 
thereby  obtaining  >>  .077  millisecond  —  a  result  which  tells  us  essentially 

nothing.  However,  Ag  without  further  experimental  or  theoretical  evidence  it 
would  be  unduly  pessimistic  to  be  concerned  about  this  possibility.8 

Proceeding  then  from  the  above-deduced  conclusion  that ,  under  atmospheric 
conditions 


>>  1.3  millisecond, 


g 


let  us  discard  the  electronic  excitation  energy  of  1L  -  0-  states  as  a  source  of 

8  ^ 

atmospheric  heating.  The  fate  of  the  vibrational  excitation  remains  to  be  con¬ 
sidered.  Here,  the  engineering  situation  is  not  favorable.  Specifically,  as 
will  now  be  shown,  it  turns  out  that  the  presence  of  water  vapor  in  concentrations 
encountered  under  typical  mid-latitude  conditions  leads  to  vibrational  de¬ 
excitation  of  02  in  times  substantially  less  than  one  millisecond. 


The  required  information  concerning  collisional  quenching  of  vibrational 
excitation  in  02  is  obtained  from  experiments  on  absorption  and  dispersion  of 
ultra-sound  in  gases.  These  experiments  yield  the  so-called  "Napier"  relaxation 
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frequency  fN  *  1/2  it  t^,  where  t^,  the  "Napier"  time,  is  the  time  required  for 
thermal  equilibration  between  vibrational  and  translational  degrees  of  freedom  of 
a  gas;  tn  may  be  identified  as  the  lifetime  of  vibrational  excitation,  and  is 
hence  of  direct  interest  to  us. 


o 

The  results  of  Henderson,  Clark,  and  Lintz  will  form  the  basis  of  the 
present  discussion.  According  to  these  authors,  tn  in  (^-^O  mixtures  is  given 
by  the  expression 


~ —  -  (7  +  18  3h  +  132h2)p 

2tttn 

(with  1/t  in  sec  ^) ,  where  p  is  the  pressure  of  the  main  constituent,  0_,  in 
^  -  3 

atmospheres,  and  where  h  -  10  N„  n/Nn  .  One  notes  that  when  h  exceeds  unity 
(N„  /N„  >  .1%),  the  term  quadratic  in  the  water  vapor  concentration  becomes 

HoO  ^2 
dominant . 

Before  going  on  to  evaluate  the  engineering  situation,  let  us  digress 
briefly  to  discuss  the  mechanism  of  the  quadratic  term.  The  currently  believed 
mechanism  (first  proposed  by  Tuesday  and  Boudart^)  consists  of  a  two-step  process. 
The  first  step  is  the  transfer  of  vibrational  energy  from  0£  to  ^0  in  a  binary 
collision;  this  transfer,  by  virtue  of  the  close  energetic  proximity  of  the 
lowest  vibrational  level  of  0„(1556  cm  to  the  lowest  bending-mode  excitation 
of  ^0  (1595  cm  ),  is  quasi-resonant,  and  hence  assumed  to  be  rapid  enough  to 
establish  a  common  vibrational  temperature  for  the  two  species,  i.e., 


9 


“ft 

where  N_  and  N„  n  denote  the  concentrations  of  vibrationally  excited  molecules 
u_  n_U 

of  O2  and  H2O.  The  second  step  consists  of  the  conversion  of  ^O-vibrational 
excitation  into  translatory  (or  rotatory)  energy  via  ^0  -  H20  collisions.  This 
reaction,  although  quite  rapid  compared  to  other  bimolecular  vibrational  de¬ 
excitation  processes  (e.g.,  0„  with  itself)  is  still  taken  to  be  slow  compared 

^  £  £ 

to  the  reverse  of  the  first  step  (i.e.,  HgO  +  (^  "*  HjO  +  C^).  It  is  then  readily 
seen  that  the  quadratic  term  in  the  above  expression  for  tn  is  given  by  the 
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formula 


h2o 


N, 


Nn  'H.O  <VH„0>  Qq* 


Where  NH20  is  the  concentration  of  normal  H20  molecules,  <vH  Q>  is  some  average 
relative  velocity  of  two  colliding  H,,0  molecules,  and  Qq  is  the  cross  section  for 
vibrational  quenching  of  H20  in  such  a  collision.  Inserting  into  this  relation 
the  experimental  value  for  1/t2,  namely  2ir(132)h2  sec"1,  and  using  the  formula 

h  "  10  NH20/N02’  one  finds  hh20  <vH20>  Qq  *  8.3  x  108  sec-1 /atmosphere. 

With  V  "  2,7  x  lol9cn1"3  at  one  atmosphere,  and  with  <vH  _>  ~105cm/sec_1,  one 
has  2  ^  .  H2° 


Qq  —  3  x  10"16cm2, 

which  is  1/10  a  gas-kinetic  cross  section.  This  is  perhans  one  order  of  magni¬ 
tude  higher  than  the  usual  rate  of  vibrational  deactivation  in  reactions  involving 
polyatomic  molecules.11 

It  should  finally  be  stated  that  the  basic  reason  for  the  credibility  of 
the  Tuesday-Boudart  mechanism  (as  compared  to,  e.g.,  a  three-body  interaction 
between  one  oxygen  and  two  water  molecules)  is  that  the  quadratic  term  is  unique 
to  H20;  D20  and  HDO  give  rise  only  to  terms  linear  in  their  respective  concen¬ 
trations.  Since  the  vibrational  frequencies  of  these  molecules  are  not  close  to 
that  of  02,  the  absence  of  the  quadratic  term  finds  a  natural  explanation  in  the 
Tuesday-Boudart  theory.12 

Having  discussed  the  mechanism  of  the  quadratic  term  (the  linear  term  being 
presumably  associated  with  the  reaction  0*  +  H20  +02  +  H20  +  K.E.,  let  us  evalu¬ 
ate  the  engineering  situation.  For  this  evaluation  it  is  useful  to  present  a 
table  of  h  as  a  function  of  altitude  for  an  average  mid-latitude  location.13 

Altitude  (km)  Mixing  Ratio (W  =  M^p/M^)  h  -  W(29/18)  x  103 


0 

.63 

X 

CM 

1 

o 

cH 

10.2 

2 

.38 

X 

10"  2 

6.1 

4 

.19 

X 

10" : 2 

3.1 

6 

.09 

X 

10-2 

1.45 

8 

.027 

X 

IQ-2 

.43 
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With  h  at  sea  level  equal  to  10.2,  it  is  immediately  apparent  that  the  quadratic 
2 

term  2tt  (132)h  is  dominant  and  yields  a  value 
1  4-1 

— ~  «*8.5  x  10  sec  >>/inverse  millisecond 


* 

for  the  deactivation  rate.  Furthermore,  the  corresponding  lifetime  t_  remains 

u2 

less  than  a  millisecond  for  altitudes  below  a  value  **6.3  km.  It  is  thus 
apparent  that,  for  engineering  purposes,  one  must  consider  the  vibrational  part 
of  the  absorbed  laser  energy  as  contributing  its  full  quota  to  atmospheric 
heating. 

The  energy  bookkeeping  now  proceeds  as  follows  5  if  v  ■  1/A  is  the  wave- 

4  4  “In 

number  of  the  laser  beam,  one  observes  (with  v  1  ■  10/1.26  *  .794  x  10  cm  ) 

61  14 

that  the  fraction  of  energy  available  for  atmospheric  heating  is 

v  ~  vel  _  v-7940  . 

v  v 

Taking  v  equal  to  the  Penner  value  (Appendix  E  of  the  above-quoted  Summer 
Report),  10^/1.047  ■  9550  cm  \  one  has 


v 


Hence,  for  the  effective  abaorption  cross  section  (i.e.,  the  cross  section  for 
effective  conversion  of  laser  energy  into  heat),  one  has 

v  -  v 

k  ..(v)  -  -----  k(v)  -  .17  k(v)  , 

eff  v 

where  k(v)  is  the  actual  absorption  cross  section.  Using  the  value  mentioned  at 
the  beginning  of  this  report,  one  has 

k  e.(v)  =  2.5  x  10_10cm_1. 

eff 

It  is  thus  clear  that  if  it  is  possible  to  reduce  kgff(v)  another  factor  of 
ten,  the  engineering  situation  will  begin  to  look  optimistic.  One  possibility 
for  such  a  reduction  will  now  be  discussed  briefly. 
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Specifically,  it  appears  that  apart  from  the  (02)2  absorption,  Penner's 
wavelength  choice,  10,468  X,  is  not  uniquely  indicated.  Preliminary  inspection 
of  the  Babcock  and  Moore  solar  atlas^  as  well  as  more  recent  data^  from  the 
Naval  Research  Laboratories,  indicates  the  essential  absence  of  atmospheric 
absorption  lines  down  to  about  10,300  ft.  Moreover,  Penner's  calculation  for 
10,468  ft  indicated  that  the  dominant  absorption  (apart  from  the  (Oj^  band  is 
due  to  the  relatively  distant,  strong  HjO  lines  between  11,103  ft  and  11,601  ft. 

This  contribution  will  certainly  not  increase  with  decreasing  wavelength.  On 
the  other  hand,  according  to  Fig.  1  of  Ref.  3,  che  contribution  of  (^2^2  a^’8orPt^on 
diminishes  substantially  with  decreasing  wavelength  over  the  interval  10,468  ft  -* 
10,300  ft.  The  possibility  of  obtaining  the  required  factor  of  ten  in  the  wave¬ 
length  region  10,300  -  10,370  ft  appears  at  the  present  time  to  merit  further 
investigation. 
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where  R  =  1  cm  is  the  tube  radius,  the  "diffusion"  mean  free  path,  and 

vth,  a  typical  thermal  velocity  of  Oj  molecules.  Taking  for  the  diffusion 

cross  section,  QD  =  1/NJ^,  a  characteristic  gas-kinetic  value  of  3xl0-15cm2, 
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I.  INTRODUCTION  AND  BACKGROUND 

One  of  the  processes  by  which  light  may  be  removed  from  a  beam  traversing 
the  atmosphere,  and  converted  directly  to  heat,  is  collision-induced  absorption. 
During  a  collision  between  two  molecules  a  transient  dipole  moment  may  be  in¬ 
duced  by  the  overlap  of  the  electron  clouds,  and,  in  some  cases  by  a  permanent 
quadrupole  moment  on  one  of  the  molecules.  Because  of  the  induced  dipole 
moment  an  absorptive  transition  can  occur  which  is  forbidden  for  the  isolated 
molecule.  The  transition  may  be  strictly  translational,  or  it  may  have  both 
a  translational  and  an  internal  (rotational,  vibrational  or  electronic)  com¬ 
ponent.  The  former  case  corresponds  to  an  instantaneous  transformation  of  all 
the  photon's  energy  into  heat.  In  the  latter  case,  if  the  photon's  energy  ex¬ 
ceeds  the  energy  needed  for  the  internal  transition,  then  the  excess  goes 
instantaneously  into  translational  motion,  and  therefore  into  heat.  Conversely, 
if  the  photon's  energy  is  less  than  that  needed  for  the  internal  transition, 
the  difference  is  provided  from  the  translational  energy  of  the  molecules; 
this  process  corresponds  to  instantaneous  cooling.  In  either  event  the  in¬ 
ternal  energy  is,  in  general,  eventually  transformed  into  heat  by  collisional 
relaxation. 

Coll is ion- induced  absorption  is  not  to  be  confused  with  ordinary  pressure- 
broadened  absorption.  Ordinary  pressure  broadening  applies  to  internal 
transitions  which  are  permitted  for  isolated  molecules;  the  collisions  result 
iii  a  modification  of  the  spectral  shape  of  the  line,  but  do  not  affect  the 
magnitude  of  the  frequency-integrated  absorption  coefficient,  which  is  pro¬ 
portional  only  to  the  density  of  molecules  of  the  species  making  the  transition. 
The  integrated  absorption  coefficient  for  collison- induced  absorptions,  how¬ 
ever,  is  proportional  to  the  product  of  densities  of  the  two  molecular  species 
involved  in  the  event  (or  to  the  square  of  the  density,  if  only  one  species  is 
involved).  At  standard  pressure  ordinary  absorption  usually  dominates  strongly 
over  collision-induced  absorption  in  the  spectral  vicinity  of  a  permitted 
transition.  However,  in  the  neighborhood  of  a  forbidden  transition  collision- 
induced  absorption  becomes  the  dominant  process. 

Collision-induced  absorption  is  characterized  spectrally  by  broad, 
continuous  bands  having  widths  (in  reciprocal  wavelengths)  of  the  ort’.r  of 


54 


-1  -2 

100  cm  (about  10  ev  In  energy  unite) .  These  widths  may  be  regarded  as 

following  from  the  uncertainty  principle  because  of  the  transient  nature  of 

the  Induced  dipole  moments,  which  exist  only  during  the  collision  period  of 
-14 

several  times  10  seconds.  The  bands  are  centered  at,  and  superposed  on 
forbidden  rotational  bands,  rotational-vibrational  bands,  etc.  In  addition 
there  are  bands  corresponding  to  purely  translational  transitions,  which  peak 

at  /’*'100  cm  ^ .  At  their  peaks  the  absorption  coefficients  range  in  magnitude 

*  _6  _i  _2  -4  -1  -2 

from  the  order  of  10  cm  amagat  to  the  order  of  10  cm  amagat 

A  number  of  collision-induced  absorption  bands  have  been  observed  ex¬ 
perimentally,^  ^  and  the  corresponding  absorption  coefficients  have  been 
measured  in  the  neighborhood  of  the  peaks.  The  mean  free  paths  are  quite  long 
at  standard  pressure,  so  the  measurements  must  be  made  at  pressures  of  about 
100  atmospheres  and  higher.  Most  measurements  are  within  about  400  cm  ^  of 
the  peak.  In  the  far  wings  none  have  yet  been  achieved. 

Our  interest  in  atmospheric  heating  via  collision-induced  absorption 
stems  from  the  possibility  of  its  leading  to  a  self-defocuslng  of  a  laser  beam 
propagating  through  the  atmosphere.  The  defocusing  is  brought  about  by  a  de¬ 
crease  in  the  local  index  of  refraction  as  a  result  of  laser  heating;  the  net 
effect  is  an  undesirable  decrease  in  the  transmitted  power  density.  Given  an 
upper  limit  on  the  acceptable  decrease  in  transmitted  power  density  one  can 
determine  a  corresponding  "critical"  value  of  the  reciprocal  mean  free  path 
for  heating,  which  will  lead  to  the  self-defocusing  instability.  This  critical 
value  has  been  estimated  elsewhere  for  a  contemplated  choice  of  laser  power, 
frequency,  and  pulse  duration. 

The  peak  values  of  reciprocal  mean  free  path  for  collision-induced 
absorption,  which  obtain  at  atmospheric  pressure  for  the  more  abundant 
constituent  gases  like  Nj  and  0^,  range  from  '•*'10  ^cm  to  *-'10  cm  ,  and  are 
much  greater  than  the  particular  critical  value  referred  to  above.  It  is  only 
in  the  far  wings  of  these  bands  that  the  absorption  decreases  to  the  critical 
value.  Therefore,  to  determine  whether  or  not  a  given  collision-inductd 
absorption  band  can  cause  a  defocusing  instability,  one  must  know  the  fir  wing 
behavior  of  the  band.  However,  the  far  wings  have  so  far  been  too  difficult  to 

An  amagat  is  a  partial  pressure  of  one  atmosphere. 
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study  experimentally  even  at  high  pleasure .  so  that  It  Is  desirable  to  ette.pt 
a  theoretical  prediction  of  their  behavior . 

Whether  or  not  collision- induced  absorption  presents  a  potential  limi¬ 
tation  to  laser  beam  transmission  through  the  atmosphere  depends  upon  the 
laser  power  being  considered.  The  maximum  absorption  coefficient  for  heating 
which  can  be  tolerated  (insofar  as  beam  self-defocusing  is  concerned)  is 
smaller  for  larger  laser  power.  For  some  beams  which  one  might  wish  to  use, 
such  as  the  kilowatt  c-w  beam  in  the  10  micron  wavelength  range,  obtainable 
with  the  C02  laser,  collision-induced  absorption  is  probably  not  a  serious 
problem.  For  megawatt  beams,  however,  the  critical  absorption  coefficient  is 
much  less  than  the  peak  values  of  typical  collision-induced  absorption  co¬ 
efficients;  in  this  case  the  relative  spectral  locations  of  the  laser  beam  and 
the  absorption  bands  are  crucial.  If  for  a  given  laser  frequency  and  intensity 
the  critical  absorption  coefficient  proves  to  be  exceeded,  then  collision- 
induced  absorption  has  associated  with  it  a  difficulty  not  present  for  most 
other  heating  mechanisms.  In  other  mechanisms  the  absorbed  lig.it  energy  is 
temporarily  stored  in  some  internal  state  of  the  absorbing  molecule,  and  is 
released  into  heat  energy  ;ver  a  period  of  some  collisional  relaxation  time. 
Therefore,  by  making  the  laser  pulse  duration  less  than  the  collisional  re¬ 
laxation  time  one  can,  at  least  in  principle,  overcome  the  defocusing  in¬ 
stability  resulting  from  most  heating  mechanisms.  However,  in  the  case  of 
collision-induced  absorption,  all  or  a  large  part  of  the  photon  energy  may  go 
instantaneously  into  heat,  and  the  heating  problem  cannot  be  surmounted  by  using 

a  short  laser  pulse. 

The  specific  list  of  laser  beams  (with  their  powers,  frequencies,  and 
pulse  lengths)  for  which  we  with  to  evaluate  the  capability  of  atmospheric 
transmission  has  not  yet  been  made  final,  and  in  fact  cannot  be  predicted  much 
better  than  the  characteristics  of  lasers  that  will  be  developed  in  the  next 
few  years  can  be  predicted.  In  view  of  this  circumstance  and  the  fact  that 
atmospheric  collision-induced  absorption  bands  are  scattered  throughout  the 
visible  spectrum  (they  presumably  occur  everywhere  that  forbidden  transition 
lines  occur),  we  have  concluded  that,  despite  the  current  emphasis  on  the 
relatively  low  power  C02  laBer,  collision-induced  absorption  bands  require 

investigation. 
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V 


Two  of  the  atmospheric  "windows,"  which  are  currently  under  consideration, 

are  located  in  the  wavelength  vicinity  of  one  micron  and  ten  microns  respec- 

3 

tively.17  The  ten  micron  window  (reciprocal  wavelength  ~10  cm  )  appears  to 
be  threatened  by  a  collision- induced  absorption  band  resulting  from  a  collision 
between  two  N0  molecules. ^2  This  band  is  superposed  on  a  rotational  band  of 
N2;  its  peak,  located  at  ^100  cm  ,  corresponds  to  an  absorption  coefficient 
of  ^2.5  x  10”6cm-1  amagat”2.  The  estimation  of  the  magnitude  of  this  ab¬ 
sorption  coefficient  at  103cm-1  has  been  made  a  specific  goal  of  our  program. 


The  ten  micron  window  is  also  endangered  by  a  vibrational -rotational 
band  of  C>2,  which  results  from  C>2  -  C>2  collisions  (and  probably  also  from 
02  -  N2  collisions).  This  band  has  a  peak  absorption  coefficient  of  3.5  x 
10”^cm-1  amagat-2,  with  the  peak  located  at  1550  cm  ^  It  is  considerably 
closer  to  the  window  than  the  N2  -  N2  bands  and  therefore  at  first  sight  might 
be  expected  to  dominate.*  However,  the  reverse  is  probably  the  case  because 
of  the  fact  that  collision-induced  absorption  bands  characteristically  de¬ 
crease  much  more  rapidly  on  the  low  frequency  side  of  their  peaks  than  on  the 
high  frequency  side.  Specifically,  the  frequency  dependence  of  the  absorption 
coefficient  on  both  sides  differs  principally  by  an  extra  factor  exp[-*(«o-u)/KT] 
on  the  low  frequency  side,  where  u)q  is  the  peak  frequency  and  id  is  the 
frequency  in  question.**  The  explanation  of  this  behavior  is  straightforward. 
The  peak  frequency  w  is  essentially  the  frequency  of  the  internal  transition. 

If  the  photon  energy  liw  is  greater  than  the  absorption  is  always 

energetically  possible,  since  the  excess  can  go  into  kinetic  energy  of  the 
colliding  molecules.  However,  if  f»u>  <  the  absorption  is  not  ^energetically 


*  Since  the  excited  vibrational-rotational  level  (which  is  located  at  about 
the  same  frequency  as  the  peak  of  the  absorption  coefficient)  is  at  a  higher 
frequency  than  the  window  (1000  cm-1)  the  initial  effect  is  a  cooling  of  the 
atmosphere  by  the  energy  equivalent  of  (1550-1000)  cm  1  -  550  cm  .  However,  if 
the  pulse  length  is  greater  than  the  collisional  relaxation  time  of  the  excited 
level,  the  entire  1550  cm-1  of  energy  will  b*?  deposited  in  time  for  the  pulse 
to  experience  a  heating  by  a  net  energy  equivalent  of  1000  cm  . 


**Thia  statement  and  the  explanation  that  follows  apply  only  to  transitions 
involving  internal  modes,  and  not  to  purely  translational  transitions. 
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possible  unless  the  kinetic  energy  E  of  the  molecules  is  great  enough  to  make 
up  the  difference,  i.e.,  unless  E  i  The  extra  Boltzmann  factor  on 

the  low  frequency  side,  then,  is  just  the  probability  that  the  colliding 
molecules  have  sufficient  kinetic  energy. 

The  far  wing  behavior  (  ~  103cnf 1)  of  the  N2  -  N2  rotational  band  is  then 
a  prime  target  of  our  program.  However,  as  a  first  step  we  have  chosen  to 
study  collision-induced  absorption  by  rare  gas  molecules  (specifically  Ne-A, 
He-A,  and  He-Ne  for  which  considerable  data  exists  *  ).  The  reason  is  that 

this  process  is  much  simpler;  there  are  no  internal  states  involved  and  the 
transition  is  purely  translational.  To  describe  a  translational  transition 
we  need  only  know  the  induced  dipole  moment  as  a  function  of  intermolecular 
separation.  However,  to  describe  a  rotational-translational  transition  we 
must  know,  in  addition,  the  dependence  of  the  induced  dipole  moment  on  the 
orientations  of  the  Individual  molecules.  Therefore  we  shall  first  study  the 
simpler  absorption  by  rare  gases,  then  use  the  rare  gas  absorption  coefficients 
in  the  far  wings  to  make  a  first,  primitive  estimate  of  the  N2  -  N2  absorption 
at  1000  cm  ,  before  we  go  on  to  attack  the  latter  more  directly. 

Considerable  theoretical  work  has  been  done  by  Van  Kranendonk  and  others 
on  the  frequency- integrated  collision-induced  absorption  coefficient.18”24 
Satisfactory  agreement  with  the  data  appears  to  be  obtained  by  using  a  model 
in  which  the  induced  dipole  moment  consists  of  two  terms.  The  first  is  a 
result  of  the  overlap  of  the  electron  clouds  of  the  two  molecules  during  the 

collision  and  of  their  consequent  distortion;  it  is  very  short  range  and  has 

-R  /  o 

an  exponential  dependence  ^e  on  the  intermolecular  separation  R,  and  is 

directed  along  the  intermolecular  axis  The  second  term  describes  a  dipole 
moment  induced  on  one  molecule  by  a  permanent  quadrupole  moment  on  the  other; 

*  ” 

There  is  no  collision-induced  absorption  by  rare  gas  molecules  of  the  same 
species,  because  the  spatial  symmetry  in  the  center  of  mass  system  precludes 
the  existence  of  an  asymmetric  quantity  like  the  induced  dipole  moment.  In 
the  case  of  diatomic  molecules  of  the  same  species,  like  N2  -  N2,  the  freedom 
of  orientation  of  each  molecule  relative  to  the  intermolecular  axis  destroys 
the  symmetry,  so  that  the  dipole  moment  can  be  induced. 
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it  Is  longer  range  than  the  overlap  term  and  goes  as  R  .  Since  rare  gas 
atoms  do  not  possess  permanent  quadrupole  moments,  we  shall  be  concerned  only 
with  the  overlap  moment. 

The  evidence  for  an  exponential  overlap  moment  also  includes  a  direct 

19 

quantum  mechanical  calculation  by  Van  Kranendonk  and  Bird,  who  find  that 
their  numerical  results  can  be  well  represented  by  an  e  dependence. 

We  are  aware  of  only  three  papers  concerned  with  the  frequency  dependent. 

25-27 

collision-induced  absorption  coefficient.  All  three  are  confined  to  rare 

25 

gas  molecules.  Levine  and  Birnbaum  describe  the  molecular  motion  by  a 

classical  straight  line  trajectory  and  take  for  the  induced  dipole  moment 

y(R)  -  y  yR  e_Y  R  (R/R) .  The  factor  yR  in  front  is  used  to  simulate  the 
o 

behavior  of  the  scattering  wave  functions:  the  molecules  are  actually 
scattered  by  a  potential  which  has  a  steep,  short  range,  repulsive  part  that 
causes  the  wave  functions  to  decrease  rapidly  for  small  R.  The  exponent  is 
chosen  quadratic  in  R,  Instead  of  linear,  because  the  former  choice  signifi¬ 
cantly  simplifies  the  analysis  and  permits  a  closed  form  expression  for  the 
absorption  coefficient.  The  parameter  yQ  appears  trivially  in  the  result  as 
a  scale  factor  multiplying  it;  the  spectral  shape  of  the  absorption  coefficient 
depends  only  on  y.  By  choosing  y  appropriately  they  are  able  to  fit  the  data 
remarkably  well  for  He-A.  For  Ne-A  the  fit  is  good  at  low  frequencies,  but 
their  theoretical  absorption  coefficient  begins  to  decrease  more  rapidly  than 
the  measured  coefficient  at  higher  frequencies.  From  their  ability  to  fit  the 
data  with  a  classical  straight  line  trajectory  they  conclude  that  the  absorption 
coefficient  is  not  sensitive  to  the  potential. 

Levine’  has  done  a  quantum  mechanical  calculation,  replacing  the  true 
scattering  wave  functions  by  plane  waves,  and  using  the  same  dipole  moment  as 
above,  with  the  factor  yR  again  having  the  purpose  of  simulating  the  true  be¬ 
havior  of  the  wave  functions  at  small  R.  The  results  are  again  in  closed  form 
and  differ  little  from  thode  of  the  classical  model  at  lower  frequencies.  At 
higher  frequencies  the  two  results  begin  to  diverge. 

97 

Tanimoto  has  done  a  strictly  quantum  mechanical  calculation  using  a 
more  realistic  model.  The  scattering  wave  functions  are  solutions,  albeit 
approximate,  to  the  Schroedinger  equation  with  an  exponential  repulsive 
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-2R/X 

potential  of  the  form  V(R)  =  e  .  The  induced  dipole  moment  is  taken 

to  be  of  the  form  recommended  by  the  work  of  Van  Kranendonk,  etc.,  namely, 

y(R)  *  (jq  e  (R/R) .  Tanimoto's  results  depend  on  three  parameters,  V  ,  X, 

and  p,  in  addition  to  the  trivial  dependence  on  y. .  He  uses  values  of  V  and 

o  o 

X  which  have  been  inferred  from  scattering  data, 28  and  chooses  p  to  fit  the 
absorption  coefficient  measurements.  Unlike  the  approaches  of  Levine  and 
Birnbaum,  Tanimoto's  analysis  requires  extensive  numerical  computations.  For 
some  special  values  of  p,  viz.,  X/p  ■  integer,  his  expressions  simplify  and 
there  is  a  radical  decrease  in  the  required  computations.  He  considers  only 
these  cases  and  chooses  the  most  satisfactory,  viz.,  p  *  X. 

It  is  difficult  to  make  a  direct  comparison  between  Tanimoto's  results  and 

those  of  Levine  and  Birnbaum  for  the  following  reason.  There  are  two  sets  of 

measurements  of  absorption  during  rare  gas  collisions.  The  earlier  measure- 

9 

ments  by  Kiss  and  Welsh,  are  in  the  near  wings,  covering  the  range  from  350 

~1  ~1  12 
cm  to  700  cm  ,  The  later  measurements,  by  Bosomworth  and  Gush,  cover  the 

peaks  of  the  bands,  and  run  from  50  cm  ^  to  400  cm  ^ .  Tanimoto's  analysis  was 

made  before  the  results  of  Bosomworth  and  Gush  were  available,  and  his 

parameters  were  chosen  to  fit  the  data  in  the  wings.  In  the  vicinity  of  the 

peaks  his  absorption  coefficients  have  reasonable  shapes,  but  are  incorrect 

quantitatively  by  as  much  as  a  factor  of  two.  Levine  and  Birnbaum  choose 

their  parameters  to  match  the  data  at  the  peaks.  They  have  pointed  out  that 

the  absorption  curves  measured  in  the  two  ranges  do  not  join  smoothly;  in 

fact  there  is  an  overlap  region  in  which  they  definitely  disagree. 

Consequently,  we  cannot  evaluate  the  accuracy  of  Tanimoto's  results  in  the 

neighborhood  of  the  peak,  or  compare  them  directly  with  those  of  Levine  and 

Birnbaum 

The  approach  we  have  chosen  is  very  nearly,  but  not  quite  identical  to 
that  of  Tanimoto.  Although  the  models  used  by  Levine  and  Birnbaum  are  much 
more  convenient,  and  actually  describe  the  experimental  observations  quite 


There  is  a  puzzling  feature  of  these  measurements.  Kiss  and  Welsh  find  an 
absorption  coefficient  in  the  wings  for  He-Ne  which  is  roughly  the  same  magni¬ 
tude  as  those  for  He-A  and  Ne-A.  Bosomworth  and  Gush,  however,  measure  the 
absorption  for  He-A  and  Ne-A  at  their  peaks  but  report  finding  no  absorption 
there  for  He-Ne. 
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well  up  to  about  400  cm  \  we  do  not  feel  confident  about  extending  their  use 
to  the  far  wings,  i.e.,  to  ~ 1000  cm-1.  First,  we  must  immediately  eliminate 
the  classical  model  from  consideration,  since  we  are  concerned  (at  reciprocal 
wavelengths  of  1000  cm-1)  with  processes  in  which  the  energy  of  the  absorbed 
light  quantum  is  about  five  times  as  great  as  the  initial  (thermal)  kinetic 
energy  of  the  colliding  molecules.  Second,  although  Levine  in  his  quantum 
mechanical  model  has  demonstrated  that  the  use  of  a  presumably  incorrect  dipole 
moment  (~>e-Y^^)  has  no  serious  consequence  at  low  frequencies,  our  goal  is 
to  extrapolate  to  higher  frequencies,  where  we  do  not  have  observations  as  a 
check.  Therefore,  we  cannot  start  with  an  a  priori  incorrect  model  and  gamble 
that  the  predictions  will  again  be  insensitive  to  it. 

Furthermore,  we  cannct  expect  the  higher  frequency  behavior  to  be  so  in¬ 
sensitive  to  the  potential,  the  latter  behavior  comes  from  the  close 
collisions,  and  in  thess  the  details  of  the  potential  are  more  important. 

Our  analysis  has  been  carried  to  the  point  where  it  is  now  being  prepared 
for  the  final  step:  a  numerical  evaluation  of  integrals  by  computer.  The 
differences  between  our  work  and  Tanimoto's  lie  in  the  evaluation  of  the 
integrals.  He  has  used  an  approximation  (it  will  be  discussed  in  Section  VI). 
whose  effect  it  is  difficult  to  assess,  so  that  we  feel  it  should  be  avoided 
in  spite  of  the  additional  computations  required.  In  addition  he  has  con¬ 
sidered  only  special  values  of  the  scale  size  p  of  the  dipole  moment,  so  as  to 
facilitate  the  computational  procedure.  We  plan  to  use  a  quite  accurate 
approximation  which  will  permit  us  to  perform  the  calculation  for  any  value 
of  p  without  seriously  affecting  the  complexity  of  the  numerical  program. 
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II.  GENERAL  EXPRESSION  FOR  DIPOLE  ABSORPTION 

We  begin  by  deriving  a  general  expression  for  the  rate  of  photon  absorption 
by  a  dipole.  The  absorption  and  emission  rates  are  given  by  the  following  two 
well  known  expressions, 


Rabs(a))  *  l  iT(f  |^,  e;i)  |  26(Ef-Ei-ftu))  , 

1 » f 


(2.1) 


Rennn(a,)  "  T  l  |T(w,eSf  |i)  1 26 (E^IUo-E,)  , 

1 » f 


(2.2) 


where  T(f|u>,e;i)  is  the  T-matrix  element  for  a  transition  between  an  initial 

state  i  and  a  final  state  f  of  the  colliding  atoms  accompanied  by  the  absorption 

of  a  photon  of  wave  vector  u/c  and  polarization  e,  and  Pi  is  the  statistical 

weight  factor  for  the  initial  states  of  energy  E^.  By  means  of  a  time  reversal 

argument  it  can  be  shown  that  R  (u>)  can  be  expressed  in  terms  of  the  T-matrix 

emm 

for  absorption  and  that,  as  a  result,  the  net  absorption  rate  R(io)  *  R  ,  -  R 

—  abs  emm 

may  be  written  in  the  form 


R(oi)  *  —■  l  (P1-Pf)|T(f|u,e;i)|26(Ef-Ei-1ia)). 
f » f 


(2.3) 


To  lowest  order  in  the  electronic  charge  the  T-matrix  element  is 


T (f  u,e;i) 


f2Trhll1/2  ,,  +|, 

l~J  ^  exp(ia,xj/c)  vcl*i>  • 


(2.4) 


where  \p^  and  \p^  are  the  initial  and  final  state  scattering  wave  functions  of  the 

colliding  atom  pair,  e  ,  m  ,  p  ,  and  x  are  the  charge,  mass,  momentum,  and 
th  J  J  J  J 

position  of  the  j  charged  particle  (electron  or  nucleus),  and  I  is  the  photon 
flux.  In  the  dipole  approximation  this  expression  reduces  to 


T(f  |u,e;i)  .  -i(M»l|1/2  ’ 


(2.5) 
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where 


“fi  *  (*£l?  • 


(2.6) 


_ 

and  i  e.x  is  the  total  dipole  operator  of  the  system, 
j  J  J 

We  assume  that  the  statistical  weight  factors  and  are  independent  of 
direction.  The  angular  integrations  In  the  sum  over  final  states  and  the  average 
over  initial  states  then  have  the  effect  of  averaging  over  directions  of  the 
polarization  vector  e.  Bearing  this  in  mind,  and  substituting  the  expression 
(2.5)  for  T(f|u>,e;i)  into  (2.3)  for  R(u>)  we  obtain 


4i  t,i  ffi-viif/'Wf-v*) 

1  9  * 


(2.7) 


Eq.  (2.7)  is  the  starting  point  for  many  of  the  theoretical  papers  on 
collision  induced  absorption.  We  shall  be  concerned  with  the  absorption  co¬ 
efficient  or  reciprocal  mean  free  path  for  photon  absorption  A(uj),  which  is 
related  to  the  absorption  rate  by 


A(uj)  ■  nn'VR(uj)/l, 


(2.8) 


where  n,n'  are  the  number  densities  of  the  two  species  of  colliding  atoms,  and 
V  is  the  quantization  volume.  The  statistics  are  described  by  the  Boltzmann 
distribution. 


’i  f  s  P(Ei  f)  *  const,  x  exp(-Ei  f/KT), 


(2.9) 


and  the  delta  function  in  (2.7)  requires  =  E^  +  tin),  so  that 


Pf  =  Pt  exp (-tlw/KT)  . 


(2.10) 
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Thus  we  find  for  A(oj)  , 


A(oj) 


^3  <*>[l  ”  exp(-WKT)] 


l  P1|wfi|26(Ef-Ei-tiU)  . 
i f 


(2.11) 


For  hoi  «  KT  the  dominant  behavior  of  the  purely  translational  bands  is 
provided  by  the  factor  ui[l  -  exp (-Tiai/KT)]  «  u>2,  and,  although  it  is  not  obvious, 
the  remaining  factor  in  A (<*>)  approaches  t.  nonzero  value  as  oj  -*-0.  Thus,  although 
the  purely  translational  bands  may  be  thought  of  as  centered  about  "the  line  at 
zero  frequency,  they  actually  have  their  peaks  at  some  nonzero  frequency  because 
of  this  factor  uj  .  The  term  exp(-fiu/KT)  may  be  traced  back  to  Eq.  (2.2),  where 
it  is  seen  to  be  a  reflection  of  the  stimulated  emission  rate  into  the  photon 
beam;  the  stimulated  emission  rate  tends  to  cancel  the  absorption  rate,  and  in 
fact  does  so  identically  at  zero  frequency.  At  photon  energies  much  greater 
than  KT,  such  as  correspond,  for  example,  to  vibrational  bands,  the  effect  of 
stimulated  emission  is  evidently  negligible. 


III.  PARTIAL  WAVE  EXPANSION 

In  the  Bor.i-Oppenheiraer  approximation  the  total  wave  function  <J>(R,r)  is 
wrl tten 


<KR,r)  -  x(R)  MR.r). 


(3.1) 


where  ^  is  the  internuclear  coordinate  and  r  represents  all  the 
electronic  coordinates,  x(^)  is  assumed  to  satisfy  a  Schroedinger  equation 
with  an  intermolecular  potential  V(S),  and  $$,r)  is  assumed  to  have  only  a 
slowly  varying  dependence  on  5.  The  dipole  matrix  element  y^  of  (2.6)  is 
then  given  by 

;fl  -  A .  <3-2) 

where 

1(5)  -  juj.jij)  U(*.?)|2  d3r  (3.3) 


is  the  expectation  value  of  the  dipole  moment  operator  for  a  fixed  internuclear 
coordinate  It.  It  is  understood  here  that  the  electronic  state  of  the  system 
(which  is  the  only  internal  state  for  rare  gas  atoms)  remains  unchanged  in  the 

transition. 

Since  we  are  concerned  only  with  rare  gas  atoms,  whose  properties  are 
independent  of  orientation,  we  assume  the  dipole  moment  y  to  be  dependent  only 
on  the  magnitude  of  the  interatomic  distance  R,  and  to  be  directed  along  the 
radius  vector  R  -  i£/R, 

y (R)  =  y (R)R  •  (3,4) 
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We  shall  use  two  center-of-mass  coordinate  systems  ,  one  with  Its  polar 
axis  along  the  Initial-state  wave  vector  k.  of  the  colliding  atoms,  and  the 

^  4 

other  with  Its  polar  axis  along  the  final-state  wave  vector  k^.  Assuming  that 
X.  and  Xr  are  scattering  states  for  a  spherically  symmetric  interaction  we  may 

*  *  -4  4 

expand  them  in  Legendre  functions  in  the  k^  and  k^  coordinate  systems,  re¬ 
spectively 


1  o  uii/^ 

xi  -  /=-  j  (24+1> 1  exp(16u)“kpr“  pJcos0i>  * 

(3.5) 

x'  '  #  i (2,+l)  1  *,tp(14«)  v  f  ’ 


where  V  is  the  quantization  volume  and  and  *«  are  the  standard  scattering 

phase  shifts.  The  coefficients  /Tv  (214*1)  i£  exp (16  /k)  have  been  chosen  so 

*  29 

that  the  wave  functions  have  the  correct  asymptotic  form, 


exp(i£*$)  4-  f(8) 


sxp (ikR) 


as  R  » 


(3.6) 


provided  the  radial  wave  functions  u . .  and  u,.  are  normalized  so  that 


u^(R)  sln(kR  4*  iir/2  4*  6^),  as  R  •+■  " 


(3.7) 


fi 


Inserting  the  partial  wave  expansions  (3.5)  into  the  expression  (3.2)  for 
,  and  using  (3.4)  for  w($),  we  obtain 


fi 


Vkfki 


l  (2tfi) (2**41) (-i)li1'  exp[i(6ljlt-6fJl)]  MU,i')  f(t,t'), 


1,1* 


(3.8) 


*  The  energy  absorbed  by  the  center  of  mass  in  the  transition  is 

(1iv/Mc2)  (•hu/2)  —lO-ll-ttw,  which  is  so  small  compared  to  its  initial 
energy  (-KT)  that  its  motion  may  be  completely  ignored. 
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where  M(£,£')  is  the  matrix  element 


MU,*')  - 


uf£(R)y(R)  ui£, (R)dR  , 


(3.9) 


and  !(£,£')  is  the  integral  over  solid  angle 


KM')  = 


P£(cos0f)  P£, (cos61)R  dQ1 


(3.10) 


The  evaluation  of  M(£,£')  is  reserved  for  the  next  section,  in  which  we  choose 
a  specific  model  for  p(R),  and  a  specific  scattering  potential  (so  that  the 
radial  wave  functions  may  be  obtained  by  solving  the  Schroedinger  equation) . 
I(£,£')  may  be  evaluated  straightforwardly  using  well  known  properties  of  Legendre 
functions.  The  procedure  is  outlined  in  the  appendix;  the  result  is 


-1 


P£(cos0fi)  (cos<f>flx  +  sin<|> 


6(£+l.£ 

(2£+l)  "  (2£+l) (2£+3) 


+  P, 


fi-  ;[_(2£-l)(2£+ 

(cos0f1)  z|_(2£-l)(2£+l)  (2£+l)(2£+3)  J  ’  (  *  ' 


where  0^,  are  the  polar  and  azimuthal  angles  of  k^  in  the  k^  coordinate 
system,  x,y,z  are  unit  vectors  in  the  same  coordinate  system,  with  z  in  the 
k^  direction,  and  6 ( £ , £  * )  is  a  Kroneker  delta. 

Now,  in  the  expression  (2.11)  for  the  absorption  coefficient  A(u)  the 
integration  over  final  states  contains  an  integration  over  the  solid  angle 


fi1 


d  n 


variable  Therefore,  we  shall  ultimately  need  the  integral 

If  we  use  (3.8)  for  pfl,  insert  (3.11)  for  I(£,£’),  and  utilize , the  well 
known  orthogonality  relations  for  Legendre  functions,  we  obtain 


fi* 


fi1 


d  Q 


fi 


l  {(£+1)  |M(£+1,£)|2+  £|M(£-1,£)|2}  .  (3.12) 


.,2,  2.  2  . 

V  k^k^ 


D/ 


As  the  quantization  volume  is  made  arbitrarily  large  the  summations  over 
states  become  integrals  according  to 


I 


(2tt) 


d  kj,  as  V  +  ®  . 


(3.13) 


Thus  the  summation  over  final  states  in  (2.11)  becomes 


mk. 


I  |y f  |2  5(E  -E-M-  — —L 
f  (2ir)J  fi2 


|ufil  d^fl  » 


(3.14) 


where  m  is  the  reduced  mass  of  the  atom  pair,  and  comes  from  the  density  of 
final  states. 


The  normalized  Boltzmann  distribution  is 
3  /  2 

Pi  “  (  V*  ^2mdCrj  eXp(-Ei/KT)  * 


(3.15) 


Combining  this  expression  with  (3.14)  and  (3.12)  in  Eq.  (2.11)  for  A(ui)  we 
finally  find,  in  the  limit  of  infinite  quantization  volume, 


A(io) 


3^  T"  (zSf)  “t1-  «*P(-'»/CT)] 

^  ^  |  |M(J',£)|2  exp(-E/KT)  dk  , 


(3.16) 


£=0  £'-£±1 
£'*0 


where 


E  -  ■ft2k2/2m,  E'  -  ft2(k,)2/2m  , 


E'  ■  E  +  liu)  . 


(3.17) 

(3.18) 
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IV.  EVALUATION  OF  THE  DIPOLE  MATRIX  ELEMENT 


We  wish  to  evaluate  the  matrix  element 


MU, I') 


uf4(R)y(R)uljll(R)dR. 


(3.9) 


To  do  so  we  must  specify  the  model  to  be  used.  First,  in  accordance  with  the 
discussion  of  Section  I  we  take  the  dipole  moment  to  be 


exp(-R/p)  R  , 


(4.1) 


so  that  y(R)  ■  y.  exp(-R/p).  Second,  we  take  for  the  potential 


V (R)  -  Vq  exp(-2R/X)  . 


(4.2) 


30  31 

Support  for  this  model  is  cited  by  Mason  and  by  Abrahamson.  More  precisely, 
the  potential  they  discuss  is  of  the  form 


fR, 


V (R)  -  Vq  exp(-2R/X)  -  Vj. 


(4.3) 


-3 

However,  Mason  gives  as  the  value  of  the  potential  at  its  minimum  6.35  x  10  ev 

_3 

for  Ne-A,  and  2.88  x  10  ev  for  He-A.  Since  we  shall  be  concerned  with  a  gas 

_2 

having  a  thermal  distribution  of  energies  with  an  average  of  —2.5  x  10  ev, 
we  may  drop  the  attractive  term  and  use  the  form  (4.2).  Therefore,  the  radial 
wave  functions  must  satisfy  the  radial  Schroedingsr  equation 


d2u. 


dR 


+  [k2  -  -  A2  exp(-2R/x]  u£  -  0  , 


(4.4) 


where  k2  *  2mE/h2,  and  A2  =  2mVo/ft2. 
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Now,  this  equation  cannot  be  solved  in  closed  form,  but  the  equation 
d2u 

£  2  2 

- 2  +  [q  -  A  exp (-2R/A) ]u  =0  . 

dRZ  1 


can  be.  We  shall  adopt  the  following  approximation,  which  is  used  by  Tanimoto,27 
and  appears  to  be  due  originally  to  Takayanagi .32  We  substitute  Eq.  (4.5)  for 
Eq.  (4.4)  with  q2  taken  to  be 


2  «  k2 

R2  ’ 


(4.6) 


where  Rq  is  the  classical  turning  point,  defined  as  the  solution  to 

k2  _  il  .  a2  exp(.2Ro/x)  .  o  .  (4.7) 

o 

The  idea  of  the  approximation  is  as  follows.  We  wish  to  evaluate  the  integral 
for  M(£,fc?)  given  by  (3.9).  The  integrand  decreases  rapidly  for  large  R  because 
of  the  expdnential  decrease  of  p(R),  so  that  it  is  not  important  that  the  ap¬ 
proximations  to  ufJl  and  ui£  be  accurate  at  large  R.  Both  the  real  and  approximate 
wave  functions  are  sinusoidal  for  large  R  (recall  Eq.  (3.7)),  but  the  details 
of  the  oscillatory  behavior  do  not  have  much  effect  on  the  integral  (3.9).  The 
integrand  decreases  rapidly  for  small  R,  because  the  wave  functions  must  turn 
down  sharply  inside  the  classically  forbidden  region  (i.e.,  for  R  <  R  ).  in 
this  region  (R  <  Rq),  as  well,  the  details  of  the  wave  functions  are  not  too 
important,  provided  the  approximations  turn  down  sharply,  as  do  the  real  wave 
functions.  The  principal  contribution  to  the  integral  then  comes  from  the 
vicinity  of  the  larger  of  the  two  turning  points,  since  the  integrand  is  washed 
out  by  the  corresponding  wave  function  for  smaller  R,  and  by  p(R)  for  larger  R. 

1hus  it  is  important  that  the  wave  functions  turn  down  at  the  right  place, 
i.e.,  that  the  turning  points  be  correct.  To  this  extent  the  approximation  is 
accurate,  since  the  turning  points  Rq  of  (4.4)  and  (4.5)  have  been  arranged  to 
be  identical.  It  is,  however,  subject  to  the  criticism  that  the  true  and  ap¬ 
proximate  wave  functions  cannot  be  guaranteed  to  turn  down  at  similar  rates 
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for  R  <  Rq;  since  p(R)  rises  rapidly  as  R  decreases,  it  is  conceivable  that  a 
wave  function  which  decreases  too  slowly  will  lead  to  an  appreciable  over¬ 
estimate  of  the  integral,  and  vice  versa. 

-R/ A 

By  making  the  substitution  z  =  AAe  one  may  convert  (A. 5)  into  the 

33 

canonical  form  of  Bessel's  equation.  The  solution,  having  the  required 
boundary  condition  u^(0)  =  0,  is 


“(<R> 


H* 


_V<A» 


l-i^ 


(A. 8) 


where  Ij.  (z)  represents  the  Bessel  function  of  imaginary  argument  (and  imaginary 

index  ip),  and  N'  is  a  normalization  constant.  From  the  asymptotic  form  of  the 
3A 

Bessel  function 


(2ttz) 


1/2 


1+I±4h1  + 

8z 


and  the  values  of  AA  and  qA  it  may  be  seen  that 


IiqA(AX)  "  X-iqA(AX)  ’ 


so  that  we  may  write 


(A. 9) 


(A. 10) 


-  n  : 


Ilql(AXe-R/A)  - 


(A. 11) 


Specifically,  from  Mason's  conclusions3^  (V  •  1.07  x  10  ^erg,  A  •  A. 86  x  10-^cm 

-9  0  —9 

for  Ne-A,  and  Vq  -  2.09  x  10  erg,  A  -  5.28  x  10  cm  for  He-A)  it  follows  that 

for  Ne-A:  AA  -  3.16  x  103  and  (recall  Eq.  (A. 6))  qA  *  kA  =  7.6  (at  thermal 

energy);  and  for  He-A:  AA  *  7.95  x  10^  and  qA  ^  kA  ■  A. 3  (at  thermal  energy). 

Therefore,  the  first  term  of  the  asymptotic  expansion  (A. 9)  is  suf-ficient,  and 

(A. 11)  follows. 
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From  the  series  expansions  of  the  Bessel  functions 

tic 

U/2 


33 


that  u^(R)  has  the  asymptotic  form 


u^(R)  +  -2 IN 


sin[qR-qXln(AX/2)  +  «*»  (qX)]  , 


it  can  easily  be  seen 


as  R  +  »  ,  (4.12) 


where  <)i  is  the  phase  of  the  gamma  function  of  complex  argument, 
r(l+iy)  -  |r(l+iy)  |ei(|,(y). 


(4.13) 


By  comparison  with  the  required  asymptotic  form  (3.7)  it  may  then  be  seen  that 
the  normalization  constant  N  must  be  chosen  to  be 


N 


_i  TrqX 
2  [sinhirqX 


(4.14) 


Finally,  noting  that  the  modified  Bessel  function  of  Imaginary  argument  (and 

35 

imaginary  index)  is  defined  by 


-IT 


Kiy^  21  sinhiry  1‘iy 


m 

9 


we  may  write  u  in  the  form 


ufc(R)  "  sinhfrqX)*^  K^^^(AXe  R^)  . 


(4.15) 


(4.16) 


Inserting  wave  functions  of  the  form  (4.16)  into  the  integral  (3.9)  for 

-R/X 

M(Jt,2,' ) ,  and  making  the  substitution  z  *  AXe  ,  we  now  obtain 


,.^1/2 


,1/2  _ i. 


M(i,£')  ■  -^-(Trq'X  sinh  7rq *  X)  (irqX  sinh  TTqX)'  w 

(AX)A/P 


F(AX) ,  (4.17) 


F(AX) 


AX 


Kiq’x(z)  z>‘h'1  dz  ’ 


(4.18) 
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where  we  have  now  modified  the  notation  so  that  the  unprimed  variables  refer  to 

the  initial  state  i  and  the  primed  variables  refer  to  the  final  state  f.  For 

36 

large  argument  the  modified  Bessel  functions  go  asJ 

K-v(z)  ~  Or/2z)1/2  e“Z  ,  (4.19) 

3 

so  that,  since  AX  —  10  ,  we  may  extend  the  limits  of  the  integral (4. 18)  to 
infinity  with  negligible  error.  The  resulting  integral  can  be  evaluated  in 
closed  form, ^ 

_X/p-3  2  o 

F(-)  -  r(A/p)  lni/2  +i(q,+q)l/2]|  |  r[A/2p+i(q'-q)A/2]  T  •  (4.20) 


V.  FINAL  EXPRESSION  FOR  THE  ABSORPTION  COEFFICIENT 

We  now  collect  the  results  of  the  previous  sections  and  write  the  final 
expression  for  the  absorption  coefficient  A(w).  It  is  given  the  following 
integral  over  the  thermal  energy  distribution  of  the  colliding  molecules,  and 
summation  over  their  initial  and  final  angular  momenta, 

00 

,  .  .  •  -WKT.  V'  V' 

A(uj)  *  Ew(l-e  )  )  [ 2 

i-O  l'-i*l 
0 


x 


|  F[X/2p+i(q'+q) X/2] |4  | r[X/2p+i(q'-q)X/2] |4 


X  exp(-Ti2k2/2mKT)  sinh  irq'X  sinh  itqX 


SSi 

k' 


dk, 


(5.1) 


where  the  constant  B  is  given  by 


2,1  . 

3/2  ,  , 

X/p 

2s  - 

fl2 

(  2*2 

1 

3  1>2c 

27rmKTX2 

[mV  X2 
'  o  J 

r2(x/p) 

the  final  state  wave  number  k*  la  a  function  of  the  Initial  state  wave  number  k, 


k' 


xl/2 


+  2mu)/ti 


(5.3) 


the  initial  and  final  atate  turning  points  R,  and  ^  must  be  determined  for  all 
values  of  the  pairs  (M>  and  <k\i'),  respectively,  as  the  solutions  to 


k 


2 


k' 


2 


2mV 

o 


exp(-2RQ/X) 

exp(-2R^/X) 


l' ( l'+l ) 

n.2 


(5.4) 
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so  that  R^  and  R^  are  ultimately  functions  of  the  integration  variable  k  and 
the  summation  variables  i  and  l' t  and  finally  q  and  q'  are  defined  as  functions 
of  k,  l,  £'  by 


k2  _  £(*+!) 


1/2 


k,2_  l'(i'+l) 


,2 


1/2 


(5.5) 
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VI.  REMARKS  ON  THE  EVALUATION  OF  THE  INTEGRAL 

Eq.  (5.1)  for  the  absorption  coefficient  A(u>)  is  equivalent  to  Tanimoto' s 
expression.27  The  differences  lie  in  approximations  used  in  numerically  evalu¬ 
ating  the  integral.  Tanimoto  has  done  machine  computations  for  reciprocal 
wavelengths  up  to  uj/2ttc  -  500  cm-1,  and  has  developed  an  analytical  approximation 
which  he  asserts  duplicates  the  machine  results  quite  closely  above  300  cm  . 
Since  we  plan  to  do  an  independent  computer  evaluation  of  A(u>) ,  which  is  likely 
to  be  nontrivial,  we  feel  an  obligation  to  indicate  why  we  do  not  instead  simply 
use  Tanimoto 's  results  to  predict  the  absorption  by  rare  gases. 

In  addition  to  the  mathematical  reasons,  it  would  be  desirable  to  re-do 

12 

the  calculations  and  choose  the  parameters  to  fit  the  more  recent  data  at  the 
(presumably)  more  easily  measurable  peaks  of  the  absorption  coefficients; 

9  ^ 

Tanimoto' s  parameters  were  chosen  to  fit  the  older  data  in  the  near  wings. 

Now  let  us  return  to  the  mathematics.  For  certain  special  cases  the  gamma 

functions  of  complex  argument  reduce  to  more  elementary  functions.  In  particular 
38 

we  have 


j  r  ( 1/ 2  +  iy)  j‘ 


Tr/cosh  Try, 


T(1  +  iy)  | 2  -  iry/sinh  Try, 


(6.1) 

(6.2) 


By  combining  these  with  the  general  relation  for  gamma  functions 


r(z  +1)  ■  zT(z) 


(6.3) 


one  may  significantly  simplify  the  integrals  in  Eq.  (5.1)  provided  X/p  has  any 
integral  value.  For  this  reason  Tanimoto  has  chosen  to  consider  only  integral 


This  statement  contains  the  implicit  assumption  that  the  more  recent  data  is 
also  the  more  reliable  (recall  that,  as  mentioned  in  Section  I,  the  two  sets 
of  measurements  are  not  strictly  consistent  with  one  another),  which  in  turn 
is  based  only  on  the  assumption  that  the  larger  the  absorption  coefficient 
the  more  easily  it  can  be  measured.  Unfortunately,  we  are  aware  of  no  solid 
basis  for  choosing  either  set. 
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values  of  A/p  and  has  found  A/p  -  1  to  provide  the  best  fit.  One  might  consider 
going  a  step  further  and  calculating  A(ui)  only  for  integral  values  of  A/p,  but 
interpolating  to  obtain  A(w)  for  nonintegral  values.  However,  one  is  still  left 
with  no  satisfactory  evaluation  of  A(u>)  for  the  range  A/p  <  1. 

We  have  chosen  the  following  alternative.  Stirling's  asymptotic  formula 
38 

for  the  gamma  function, 


r(z)  ~  /u  ZZ_1/2  e"Z(l  +  ^  +  ...  )  , 


(6.4) 


is  very  accurate  for  |z|  5  2.  Therefore,  by  using  (6.3)  we  may  write 


r(z) 


rU+2) 

z(z+l) 


(6.5) 


and  then  use  Stirling's  formula  in  (6.5)  for  T(z+2) .  The  result  is  (z  -  x+iy) 


|r(z) 


2  ir  1 

-(x+2)2 

x+3/2 

(*  +y  ) 

(*+i) 

2  +  A 

1  + 


x  +  2 


2  2T 
6[(x+2r  +  y 


f 


X  exp  I -2  x  +  2  +  y  arctan 


x+2 


]• 


(6.6) 


which  is  quite  accurate  for  x  >  0.  Thus  we  can  evaluate  A(ai)  for  all  positive 
values  (and  therefore  all  physical  values)  of  A/p. 

The  other  approximation  to  be  considered  has  to  do  with  the  turning  points 
R  aud  r' ,  which  are  defined  as  the  solutions  of  Eq.  (5.4)  and  therefore  are 
functions  of  k,  £,  £'.  Tanimoto  has  chosen  to  ignore  the  k,  £,  £'  dependence 
and  use  one  fixed  value  of  Rq  (and  the  same  value  for  R^) ,  namely  the  turning 
point  for  thermal  energy  and  zero  angular  momentum,  i.e.,  the  solution  to 

2  2mV 

k  thermal  ‘  ~f  e*P<-2R0/X).  (6.7) 

If  we  assume  that  the  principal  contribution  to  the  integral  in  (5.1)  does  in 
fact  come  at  the  thermal  value  of  k,  there  appears  to  still  be  a  serious 
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problem  from  ignoring  the  angular  momentum  dependence  of  Rq,  and  R'.  The  turning 

point  in  the  absence  of  the  repulsive  angular  momentum  potential  is  necessarily 

smaller  than  in  its  presence,  so  that  when  the  smaller  values  of  R  and  R'  are 

o  o 

used  in  (5.5)  the  repulsive  effect  of  the  angular  momentum  terms  is  magnified 
and  the  effective  wave  numbers  q  and  q'  are  diminished.  As  a  result  the  wave 
function  does  not  penetrate  as  far  in  and  does  not  "see"  as  much  of  the  dipole 
moment.  Thus  it  would  appear  that  this  approximation  systematically  under¬ 
estimates  the  dipole  matrix  element  M(£,£*)  of  (3.9),  and  therefore  underestimates 
A(oj)  .  It  is  true  that  for  small  values  of  £  the  angular  momentum  potential  is 
negligible  at  the  turning  point  compared  to  the  exponential  potential.  However, 
we  have  made  some  numerical  estimates  of  the  terms  contributing  to  the  summation 
over  £  in  (5.1)  for  the  special  case  in  which  the  Boltzmann  distribution  over 
energy  is  replaced  by  a  delta  function  peaked  at  thermal  energy;  we  founa  that 
the  maximum  contribution  occurred  for  £  *•»  30,  and  significant  contributions  still 
occurred  for  £  —  60. 

Thus  we  feel  that  the  approximation  of  a  constant  turning  point  is  un¬ 
satisfactory.  The  remainder  of  our  program  includes  a  solution  of  (5.4)  to  obtain 
Rq  and  R^  as  functions  of  k,  £,  and  £f . 


APPENDIX 


In  this  appendix  we  outline  the  procedure  for  evaluating  the  integral 
!(£,«,')  =  P£(cos0f)  P£f(cos6i)  R  d^  .  0 


First,  in  the  coordinate  system  with  polar  axis  along  k^,  R  may  be  written 


R  -  sin8i  cosi()i  x  +  sin01  sin^  y  +  cos61  z  . 


Next,  P£(cos0f)  may  be  written  in  terms  of  the  variables  of  this  same  coordinate 
system  by  using  the  addition  theorem  for  Legendre  polynomials, 


P£(cos0f)  -  P£(cos0fi)  P£(cos0i) 


+  2  l  P£“(cos0fi)  P£m(cosfii)  cos  m(^-4.fi),  (A3) 


-T 

where  0^  and  <j>^  are  the  coordinates  of  k^  in  the  system.  Using  the  addition 
theorem  in  (Al)  we  obtain 


■  £(£+1)  P£ 


‘U  -  2 

P£1(cos0fi)  cos«t>t  ±  P£(coe0)  P£,(cos0)  sin  0  d0. 


Iy(^»^')  *  £(£+1)  P£ 


fTT  -  2 

P£1(coS0fi)  sin^fi  ]  P£(cos0)  P£,(cos0)  sin  0  d0, 


fTT 

Iz(£,£')  =  2*  P£(cos0fi)  P£(cos0)  P£,  (cos8)  cos0  sine  d0  .  (A6) 

•  n 
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The  integrals  in  (A4,  5,  6)  may  now  be  evaluated  using  the  identities 


40 


(2n+l)  (l-z2)1^  Pn"(z)  -  p£}<«>  -  £}<*>  , 


(2n+l)  z  Pnm(z)  “  (n-m+1)  P®+1(*)  +  (n+m)  p“_j.U)  » 


and  the  orthogonality  relation  for  Legendre  functions, 
rl 


40 


*"<■>  P«  (z)dz  -  ^  ^frtt 


-1 


(£+  m 

ill 

(£“  tn 

|)T 

t  6(4,1’)  , 


where  6 (£,£')  is  a  Kroneker  delta.  The  result  is 
p/Ccose)  P0,(cos8)  sin20d0  -  (2£-l|(21+l) 


2£(£+l) -  g  /  £ ' ) 

(2£+3)(2£+l)  614+1,4  1  ’ 


TT  21  *  , .  ,  ,t\ 

P^(cos0)  P^, (CO80)  CO80  sin0d0  -  (2£_1) (2£+l)  5 ^ 4-1,4  ' 


4.  2  (£4-1)  g  /  £ •  \ 

(2£+3)(2£+l)  (  *  '  ’ 


(A7) 

(A8) 


(A9) 


(A10) 


(All) 


When 


these  are  substituted  into  (A4,  5,  6)  we  obtain  Eq.  (3.11)  of  the  text. 
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I.  INTRODUCTION 

In  this  chapter  we  shall  investigate  the  effects  of  Raman  scattering  on 
the  propagation  of  high  intensity  laser  light  through  the  atmosphere.  The 
general  problem  of  laser  light  propagation  is  an  extremely  complex  one  which 
involves  many  different  physical  processes.  In  order  to  make  any  headway  in 
understanding  the  effects  of  Raman  scattering  we  shall  consider  a  situation  in 
which  the  effects  of  self-focusing,  etc.,  have  produced  a  stable  beam  (or  fila¬ 
mentary  "hot  spot"  within  a  beam)  of  uniform  intensity  and  cross-section.  Ex¬ 
tending  the  work  of  Helwarth,1’2  it  is  possible  to  derive  a  set  of  rate  equations 
which  relate  the  intensities  of  the  various  orders  of  forward  amd  back  amplified 
Stokes  waves  to  the  laser-light  intensity.  This  set  of  coupled,  nonlinear, 
partial  differential  equations  with  varying  coefficients  appears,  at  first  glance, 
to  be  quite  intractable.  However,  by  first  examining  a  highly  simplified  version 
of  these  equations  and  then  treating,  in  turn,  each  of  the  various  complications 
of  the  simplified  equations  which  comprise  the  original  equations,  it  is  possible 
to  obtain  a  coherent  description  of  the  physical  phenomena  involved. 

The  rate  of  creation  of  Raman-scattered  photons  in  a  laser  beam  is  com¬ 
posed  of  a  term  proportional  to  the  laser  beam  intensity  (spontaneous  emission) 
and  a  term  proportional  to  the  product  of  the  laser  beam  intensity  and  the  in¬ 
tensity  of  the  Raman  photons  already  present  (stimulated  emission) .  If  the  rate 
of  creation  of  Raman-scattered  photons  parallel  to  the  beam  (i.e.,  either  forward 
or  backscattered)  is  much  larger  than  the  rate  of  destruction  of  these  photons 
due  to  Rayleigh  scattering,  then  the  intensity  of  these  Raman  photons  will  be 
amplified  at  the  expense  of  the  primary  laser  photons  until  the  primary  beam  is 
totally  depleted.  The  specific  manner  in  which  the  primary  beam  is  extinguished 
will  be  the  subject  of  sections  II  through  V.  Section  VI  will  relate  the  current 
status  of  the  author's  researches  into  certain  transient  phenomena. 
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II.  NEGLIGIBLE  BACKSCATTER 

Let  us  consider  the  propagation  of  a  beam  of  laser  light  through  a  gas  in 
which  the  amplification  of  the  backscattered  Stokes  waves  may  be  neglected  (the 
criteria  for  the  existence  of  such  a  gas  will  be  discussed  in  section  V.). 
Assuming,  to  begin  with,  that  only  the  first  Stokes  wave  is  of  any  consequence, 
the  densities  of  the  laser  photons  p^(x,t)  and  amplified  forward  scattered  Raman 
photons  (x,t)  in  the  beam  are  related  by  the  rate  equations 


JL  ,  a_ 

3t  +  C  3x 


'  P£(apr+eo)  '  Yp*  * 


(1) 


+  pJ^apr+e)  _  YPr  ' 


(2) 


where 


c3 

— cp  o(0),  B 
2vZAv  T  ° 


CpT  °Ram.  *  Y 


CpT  °Ray . 


The  quantities  and  are  the  Rayleigh  and  Raman  total  scattering  cross- 

sections  respectively,  pT  is  the  density  of  tne  gas,  Av  is  the  line  width  of  the 
amplified  wave,  and  o(0)  is  the  Raman  forward  scattering  cross-section.  The 
quantity  3  is  given  by 


8 


«T  °<0)  —  , 


(3) 


where  is  the  effective  cone  angle  for  Raman  amplification  —  8  '  will  be  dis¬ 
cussed  in  greater  detail  later  in  this  section.  The  boundary  value  conditions 
for  Eqs .  (1)  and  (2)  are 


p£(0,t)  -  Pe(t)  ,  pQ(t)  -  0  for  t  <  0,  and  pr(0,t)  -  0  .  (4) 


We  now  make  the  substitutions: 


«■  x  3  ,  3 

C  •  y'X  -  H  +  c3^ 


in  Eqs.  (1)  and  (?).  Equation  (1)  may  now  be  rewritten  as 

9P„ 


3y 


”P £ Capr+0)  ”  (y+3c-3)  P ^ 


For  nitrogen- 


(5) 


o  *  10  ^cm^  and  0  tf  jq 

Kdni*  Ray  • 

and  hence  y  >>  3O~0,  and  we  can  neglect  8q-8  with  respect  to  y  in  Eq.  (5)  and 
make  only  a  minimal  error.  The  solution  of  Eqs.  (1)  and  (2)  with  this  approxi¬ 
mation  is 


P£(x,t) 

Pr(x,t) 

where  r(x)  -  ^ 
J(x)  - 


e-r(x) 

■  p0(8)  lMOO  • 

-  PQ(s)  e_r(x)  -p((*,t) 


■X 

dx'  y(x')  and 
o 


(6) 


+  p0(s)a(x") 


The  x-dependance  of  a,  6,  and  y  is  due  to  the  altitude  dependence  of  pt  and  Av 
for  a  laser  beam  which  is  aimed  nonhorizontally  in  the  atmosphere.  ...  From  Eq.  (2) 
we  may  note  that  no  amplification  of  pf  can  occur  unless  o  pQ  >  y.  For  Nj  at 
atmospheric  pressure  this  is  equivalent  to  having  the  laser  intensity  at 
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X  -  o  (I  )  satisfy  the  inequality 
o 

'i 

I  >  60  watts/cm<‘. 
o 

We  shall  therefore  assume  the  above  inequality  to  be  valid  for  the  remainder  of 
this  work. 

Consider  first  a  horizontally  aimed  laser  beam.  For  distances  of  a  few 
tens  of  kilometers  T(x)  is  completely  negligible.  Assuming  that  pQ(t)  -  PQ 
is  a  constant  for  t  >  o  and  defining  the  dimensionless  quantities 


y 


f  V  Fr 


a 

B  Pr 


Equation  (6)  may  now  be  rewritten  as 


y(y+1) 


y  +  e 


(y+l)z 


Fr  -  y 


e(y+l)z  _  1 


y  +  e 


(y+1)  z 


(8) 


These  solutions  indicate  that  for  y  »  1  the  laser  light  will  penetrate  the  gas 
to  a  depth  of  the  order  of  zc  -  y"1in(y)  in  such  a  manner  as  to  be  essentially 
undiminished  in  intensity  (i.e. ,  F£  *  y  for  z  *  z£) .  Beyond  this  depth  the 
laser  light  intensity  will  rapidly  decrease  to  zero  and  the  intensity  of  the 
first  Stokes  wave  will  rapidly  rise  to  the  level  of  the  original  laser  beam. 

The  effective  cone  angle  for  Raman  amplification  (fiQ)  is  a  vaguely  defined 
quantity  as  it  represents  an  attempt  to  solve  a  basically  three-dimensional 
problem  by  geometrically  correcting  one  of  the  coefficients  in  the  solution  of 
the  corresponding  one-dimensional  problem.  However,  since  -  as  we  shall 


demonstrate  below  -  the  laser  Intensity  drops  from  Iq  ■  p  hvc  to  zero  is  a 

very  small  distance,  and  the  initial  intensity  I  which  will  cause  the  intensity 

3  c 

to  drop  to  zero  at  a  distance  l  ■  8  x  10  km  (the  scale  height  of  the  atmosphere) 

is  extremely  insensitive  to  our  choice  of  8  ,  the  results  we  obtain  with  the  aid 

o 

of  this  correction  factor  should  be  quite  good. 

For  a  beam  diameter  D  and  a  z£  corresponding  to  l ,  is  given  by 


For  N2  at  one  atmosphere  Av  is  the  pressure  broadened  line  width  given  by^ 


Av 

c 


.09  cm 


-1 


Using  0(0)  •  (4i0  ^  ,  v  ■  5  x  10^  sec  ^  and  taking  D  ■  10  cm  we  may 

define  a  quantity  I  which  is  that  value  of  I.(o)  =  I  for  which  I#(.i)  is  one 

C  JC  O  )C 

half  of  I  ;  1, 

o’  ’  ’ 

V°>  =  Zo  •  'c  “d  V*)  -  7'c  •  <9> 

This  quantity  may  be  obtained  from 


p  z  ■  ap  — 
c  0  c 


or  I 


■  3.90  Cc  kilowatts/cm  , 


There  i8  the  solution  of  the  equation 


e_ 
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z  ipTo(o)  fll 


6.4  x  10 


13 


(12) 
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The  solution  of  this  equation  is  £  ■  35.4  and  hence 

c 


I 

c 


140  kw/cm 


From  Eq.  (12)  we  obtain  that 


_1_  60o 
35  n 

o 


verifying  our  assertion  that  Ic  is  insensitive  to  our  choice  of  Examining 

Eq.  (8)  we  find  that  for  IQ  -  Ic,  1^  drops  from  ( . 75) Ic  to  (.25)Ic  in  a  distance 
Ax  given  by 


-  .062  or  Ax  -  500  meters. 


and  hence  the  transition  region  is  quite  narrow  and  nQ  may  be  defined  in  a 

reasonably  unambiguous  manner.  In  fact,  6 / f?o  4  2  Ax/ 1  and  hence  the  ambiguity 

in  0  results  in  a  fractional  error 
0 


1  2Ax 
35  l  " 


4  x  10 


A  plot  of  the  fractional  transmission  to  a  depth  l  versus  the  initial 
intensity,  or 

I  ,  .  1  ±  u 

»  1  y  +  <>+!>* 

versus  y  for  fixed  z,  would  result  in  a  curve  much  like  that  for  v8,  *• 
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That  is,  I  ^1  (l)  would  be  roughly  constant  and  equal  to  one  for  0  <  I  <  I  , 

0£,  *  OC 

equal  to  1/2  for  J  ■  Ic>  and  drop  sharply  to  zero  for  Iq  >  Ic>  A  plot  of  the 

total  transmission  (1^)  versus  Iq  would  therefore  yield  a  curve  which  would 

show  a  roughly  linear  rise  from  zero  to  some  maximum  value  I  ■  I  ,  and  then 

me 

fall  rapidly  to  zero.  The  maximum  transmission  to  a  depth  l  is  readily  obtained 
from  Eq.  (8)  by  setting  the  derivative  of  F^  with  respect  to  y  equal  to  zero. 

The  result  is 


m 


1  + 

m 


—  ,  I  -  3.90  i  kw/cm  , 
-io  m 


(13) 


where  £  is  the  solution  of 
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2.  ^m 

(1  -  “)e 
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i 
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(14) 


For  the  parameters  of  Interest  here 


I  ■  120  kw/cm^  and  I  ■  124  kw/cm^ 

m  o 
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"'IT, ,  — - » -  — 

Let  ub  now  consider 


with 


.  -x/i 

PT(x)  -  PT(°>e 


e  the  maximum  intensity  which  can 

Sine.  «•  C“C'med  “"ftmo.ph.r.  ».  —  that  th.  «-> 

be  transmitted  throng  the  top  beam.  We  may  therefore 

lln.  1.  virtually  unchang.d  throughout  —t  «* 

Z~  th.  lln.  width  of  th.  «UH.d  w.v.  to 


Av  -  Avp(o) 


and  b  " 


Av  (o) 
P 


broadened  Raman  lln.  width  and  tvt  1.  th.  Un« 
wh«r.  %«»  1*  the  Pt“*Ure  m  „.lng  E„.  (15)  and  <1«  »  obtain  u,  6. 
width  of  th.  original  laaar  •  t„uU.  ln  Eqa.  (6)  “d  (7)  w. 

„d  v  a.  a  function  of  x  tc.tt.rlng)  1.  n.gllgl--  A  *°°d 

fl”d  th>t  !“ /<y) « X  -  defied  for  altitude,  gr.at.r  than  on. 

approximation  to  Jty)  w  " 
kilometer .  The  result  is 

yz0f(y)  ° 

J(y)  -  7e 


where  zq  -  ^/c  and 


K1) 


M  in  Eq.  (6)  we  obtain 
Inserting  Eq.  (18)  q 


F^y)  - 


y  +  e 


yzQf(y) 
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Choosing  (Av./c)  -  1.2  x  lO-3  cm"1  (which  is  typical  for  the  line  width  for  a 

x,  5 

Ruby  Giant  Pulse  Laser  )  we  find  that 


n 


_l_ 

75  • 


For  y  «  yQ  -  -  i  in  n  -  35  km,  it  follows  that  zQf  *  zQy/i  -  z,  and  Eq.  (19) 
agrees  with  Eq.  (8)  to  within  the  approximations  made  here. 

A  calculation  of  the  maximum  transmission  of  laser  light  to  an  altitude 
y  may  be  made  in  a  manner  similar  to  that  used  for  the  analogous  quantity  in 
the  preceding  section.  The  result  is 


I  - 

m 


I 

o 


1  +  (S  -2) 
m 


i 

o 


3.90 


kw/cm‘ 


(20) 


where  £  is  solution  of 
in 


m 


(21) 


The  last  equality  is  valid  since  for  y  >  yQ  the  effective  cone  angle  for  Raman 
amplification  (no)  is  proportional  to  y*.  For  y  >>  35  km  (i.e.,  for  above  the 
atmosphere),  the  maximum  User  intensity  which  can  be  transmitted,  and  the 
laser  intensity  necessary  at  the  source  to  achieve  thut  maximum  are 


2  2 
1-27  kw/cm  ,1-28  kw/cm  • 
m  o 

It  is  of  interest  to  note  that  the  atmosphere  will  act  as  an  amplifier 
for  Raman  waves  up  to  an  altitude  of  approximately  yQ  (about  four  scale  heights) . 
This  fact  is  most  easily  understood  if  we  note  that  the  gain  is  proportional 
to  p^/Av.  It  then  becomes  clear  that  the  amplification  of  Raman  waves  up  to 
y  occurs  because  although  the  density  of  N2  decreases  with  altitude,  the  width 
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of  the  pressure  broadened  Raman  line  decreases  in  a  proportionate  manner.  The 
net  effect  is  to  produce  a  constant  gain  for  y  <  yQ  (i.e.,  as  long  as  Av  is 
dominated  by  Av  )  and  a  consequent  decrease  in  the  maximum  transmittable  power. 

In  order  to  minimize  this  effect  the  laser  line  should  be  made  as  broad  as 
possible.  If  the  laser  line  is  significantly  broader  than  the  pressure  broadened 
Raman  line,  then  o(x)  is  proportional  to  the  local  density  and  the  solution  of 
Eqs.  (1)  and  (2)  may  be  obtained  from  Eqs.  (8)  by  making  the  substitution 


Z  -*■  T 


z 

o 


(1  -  e~x/l) 


(22) 


where  t  is  the  optical  depth.  Consequently,  the  intensity  of  laser  light  trans¬ 
mitted  vertically  to  a  point  far  above  the  atmosphere  will  be  the  same  as  that 
transmitted  a  distance  of  one  scale  height  through  air  held  at  the  constant 
pressure  of  one  atmosphere  throughout  the  beam. 

Of  further  interest  is  the  fact  that  Rayleigh  scattering  is  virtually 
always  inconsequential  for  the  intensities  of  interest  here.  Its  major  effect 
is  to  provide  us  with  a  good  reason  for  changing  8q  to  6  in  Eq.  (1). 
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IV. 


MULTIPLE  STOKES  LINES 


Let  us  now  consider  a  case  in  which  the  laser  beam  intensity  is  much 
larger  than  that  necessary  to  induce  its  conversion  to  the  first  Stokes  wave 
within  a  certain  prescribed  distance.  This  first  Stokeu  wave  will  then  propa¬ 
gate  a  certain  distance  until  it  is  converted  into  the  second  Stokes  wave,  etc. 
Assuming  that  the  pressure  is  constant  along  the  beam  path  and  that  the  effects 
of  back  scattering  may  be  neglected,  there  will,  in  general  be  N  Stokes  waves 

(of  frequency  v  -  v  ,  v  -  2v . .  v.  -  Nv  ,  for  a  laser  frequency  v  and 

a  Raman  frequency  shift  v^)  which  will  attain  a  significant  Intensity  within 
a  given  depth  in  the  gas.  In  the  dimensionless  units  of  Section  II  the  photon 
densities  of  these  Stokes  waves  are  coupled  by  the  set  of  equations: 


_d_ 

dz 


-  Ft(Fl  +  1) 


_d_ 

dz 


fa(f1  +  l)  -  f1(f2  +  1) 


dz  rN-l 


i 

i 

i 
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N-2 


(FN-1  + 


fn-i(fn 


+ 1) 


_d_ 

dz 


N 


fn-i(fn 


+ 1) 


J 


(23) 


These  equations  have  been  solved  exactly  for  the  case  of  N  ■  2.  The  solution 
is  too  complicated  mathematically  for  the  reader  to  obtain  much  use  from  its 
display  here.  Basically,  it  describes  a  situation  in  which  the  laser  beam 
penetrates  the  gas  to  a  depth  z^  ■  p  In  p  in  such  a  manner  as  to  be  essentially 
undiminished  in  Intensity.  Beyond  this  depth  the  laser  light  intensity  will 
rapidly  decrease  to  zero  and  the  first  Stokes  wave  will  rapidly  rise  to  the 
level  of  the  original  laser  beam.  Between  z^  and  ■  2y  ^  £n  p,  F^  is  roughly 
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equal  to  the  constant  At  *2  the  first  Stokes  wave  suddenly  falls  to  zero 
the  second  Stokes  wave  (F^)  rises  to  v .  A  close  examination  of  Eq.  (22) 
taking  due  note  of  the  fact  that 


N 

Fp(z)  +  Z  F  (z)  -  u 

*  ,  m 


(obtained  by  adding  Eq.  (23)  and  integrating  the  result),  indicates  that  this 
process  will  continue  for  a  very  large  number  of  orders  if  v  »  N.  Hence  if 
we  wish  to  transmit  energy  in  a  beam  vertically  through  the  atmosphere,  it 
would  be  advisable  to  choose  ^  as  the  highest  frequency  of  a  window  in  the 
atmospheric  transmissivity.  If  the  width  of  the  window  is  Nv  ,  and  if  the 
intensity  of  the  laser  beam  is  the  maximum  permissible,  then  a  light  beam 
of  frequency  -  NVr  will  emerge  from  the  top  of  the  atmosphere.  On  the 
basis  of  Eq.  (22)  we  may  estimate  this  maximum  intensity  to  be  the  critical 
intensity  for  a  distance  Jl/N.  That  is 


Im  ■  3.90  (NCc)  kw/cm2 

where 


(24) 


For  a  transmission  window  of  100  X  and  a  Raman  wavelength  shift6  of  12  X  we 
have  N  -  8  and 


I 

m 


1.0 


megawatt /cm 


2 


9 
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V. 


THE  EFFECTS  OF  BACKSCATTERING 


For  the  case  of  comparable  Raman  forward  and  back  scattering  cross 

sections  in  an  Infinite  gas  of  uniform  density  along  the  ray  path  let  us  assume 

that  only  the  first  Stokes  wave  is  of  any  consequence  (this  is  true  up  to  the 

second  critical  depth  which  is  beyond  the  distances  of  interest  in  this  section) . 

The  rate  equations  relating  p  to  the  forward  (p,)  and  back  scattered  (p,  ) 

x*  £  b 

photon  densities  are 


(25) 


(26) 


°n'Vb  +  V 


(27) 


The  transient  response  of  the  above  equations  is  extremely  complex.  Let  us 
reserve  our  discussion  of  this  response  for  Section  VI  and  consider  only  the 
steady  state  solutions  (i.e.,  where  the  time  derivatives  in  the  above  equations 
vanish)  here. 

Adding  Eqs.  (25),  (26),  and  (27),  and  integrating  the  result  we  obtain 


p^(x)  +  Pf(x)  -  Pb(x)  -  constant 


-  p0  "  Pb(°)  “  P £ (L) 


(23) 
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or, 


Pb(o)  +  pf(L)  -  pQ 


(29) 


In  the  above  we  have  taken  L  to  be  a  point  just  before  the  second  Stokes  wave 
becomes  important.  Dividing  Eq.  (26)  by  Eq.  (27)  we  readily  obtain 


The  substitution  of  Eq.  (31)  into  Eq.  (29)  yields 


(30) 


(31) 


[■ 


Pb(o) 


(32) 


If  Avf  -  Avb,  then  (af/Bf)  -  (0^/6^  ,  m  a(o)/a(ir)  and  we  may  define 

a  quantity  M  such  that 


M  -  l+fpb(o)  . 

M  satisfies  the  equation 

jjO(o)/o<u)  +  M  .  2  +  „  ,  (33) 


97 


If  o(o)  =  a ( tt)  then 


pb(°)  "  pf(L)  -  ip0 

and  the  energy  in  the  original  laser  beam  is  shared  equally  between  the  forward 

2 

and  back  scattered  waves.  However,  for  I  =  100  kw/cm  the  solution  of  Eq.  (33) 
indicates  that  for 

p,  (o)  ■  .05  p  and  p,(L)  ■  .  95  p  , 

d  or  o 


That  is,  forward  scattering  will  dominate  back  scattering  as  a  mechanism  for 
extinguishing  the  laser  beam  if  the  asymmetry  in  the  cross-section  is  as  small 
as  8 7..  Conversely,  the  back  scattering  will  dominate  the  forward  scattering 
if  the  cross  section  of  the  former  exceeds  that  of  the  latter  by  as  little  as 
8%.  This  second  case  is  particularly  distressing  since  the  atmosphere  will 
then  essentially  act  as  a  mirror  -  reflecting  all  the  energy  back  at  the 
source  -  if  the  source  intensity  exceeds  some  critical  value  for  a  given  path 
length. 

If  a(o)  ■  a(ir)  then  Eqs.  (25),  (26)  and  (27)  are  easily  solved.  The 
results  are 
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where 


2 

The  critical  intensity  defined  by  Eq.  (9)  is  130  kw/cm  which  is  virtually 

2 

unchanged  from  the  case  of  purs  forward  scattering  where  I  ■  140  kw/cm  . 
However,  the  shape  of  the  I.(x)  vs.  x  and  1.(4)  vs.  I  curves  are  somewhat 

K  Xj  O 

different  for  the  two  cases:  The  distance  in  which  I  (x)  drops  from  .75  I 

to  .25  I  is  Ax  -  850  meters  versus  500  meters  for  the  case  of  no  back  scatter. 

Similarly,  for  Iq  -  Ic:  1^(4)  *  65  kw/cro^,  and  1^(4)  -  1^(4)  =  2  x  10  ^  W/cm^; 

whereas  for  the  case  of  negligible  back  scatter  if  I  ■  I  then 

2  o  c 

1^(4)  -  1^(4)  ■  70  kw/cm  .  Consequently,  the  effect  of  backscatter  is  to  soften 
the  abrupt  transitions  in  intensity  which  were  characteristic  of  the  1^  curves. 

It  is  of  interest  to  note  that  if  o(o)  ■  o(tt)  then  the  maximum  achievable 
intensity  at  a  depth  4,  and  the  source  intensity  to  achieve  this  intensity  are 

I  ■  118  kw/cm^  for  I  -  243  kw/cm^ 
m  o 

2 

versus  120  and  124  kw/cm  for  the  case  of  negligible  back  scatter.  Hence,  I  is 

essentially  insensitive  to  the  asymmetry  in  the  cross  section  but,  since 

(.5)1  is  reflected  back  into  the  source  for  o(o)  "  o(vr),  I  must  be  roughly 
o  o 

twice  as  large  to  achieve  the  same  I  in  the  symmetric  case. 

m 
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VI .  TIME  DEPENDANCE  OF  BACKSCATTERED  WAVES 

We  were  able  to  obtain  and  as  a  function  of  space  and  time 

for  the  case  of  inconsequential  backscatter  [see  Eqs.  (6)  and  (7)1.  Unfortunately, 
we  are  unable  to  repeat  this  feat  by  solving  Eqs.  (25),  (26),  and  (27).  It  is 
possible,  however,  to  solve  these  equations  for  the  case  of  inconsequential 
forward  scattering  and  draw  some  inferences  about  the  general  case  as  well  as 
obtain  some  results  which  are  not  directly  related  to  the  problem  of  atmospheric 
propagating  but  are  nonetheless  of  physical  interest. 

Neglecting  3f  in  Eqs.  (26)  and  (27),  we  find  that  pf  -  o.  The  change 
of  variable 

•  ■  t  -  *  ,  t  -  t  +  *  , 


R.(s,t)  -  p  C*,t)  ,  ^(s.t)  -  pb(x,t) 


enables  us  to  write  Eqs.  (25)  and  (27)  as 


3R„ 


j?  -  -  W  +  «  * 


(35) 


^b 

3s 


V«Rb  +  6) 


(36) 


The  general  solution  of  these  equations  is  given  by 


R£(s,t) 


UjlL 


G(t)  -  |  F(s) 


(37) 


Rjj  (s ,  t)  - 


-i  +  i 


MiL 


°  “  G(T)  -  f  F(s) 


(38) 
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where  F  and  G  are  arbitrary  functions  and  F  and  G  are  the  derivatives  of  these 
functions  with  respect  to  their  arguments. 

Let  us  now  consider  the  solution  of  Eqs.  (35)  and  (36)  subject  to  the 
boundary  value  conditions 

pJt(o,t)  -  PQ(t)  ,  p£  -  pb  -  o  for  t  <  o  .  (39) 

Equation  (36)  indicates  that  is  a  continuous  function  of  s  and  hence,  since 
Pb  ■  o  for  t  <  0  we  have 


V°>  x)  “  °  * 

Using  Eqs.  (39)  and  (40)  in  Eqs.  (37)  and  (38)  we  obtain 

-  8  £ 

„  ,  ,  .  e _ -  M(s,t) 

R£(s,t)  -  Pq(s)  1  _ 


n  .  8  M(s.t) 

V8,t)  "  a  1  -  M(s,t)  » 


where  2x  ■  c(t-s)  and 


M(s,t)  - 


dC  pQ(0  exp  (5-^) 


§  dV  P n(V) 

2  o 


.  (43) 


If  p  (t)  ■  constant  ■  p  then  for  s  >  o 
o  o 


*b  "  po 


1  -  e 


[(y+l)e  -  y]  +  ye 
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(45) 


R*  ‘ 


where  2g  ■  8 


V°> 


p  1  +  e~gs 

°  [(y+l)eZ  -  y]  +  e“gs 

+  opQ.  Equation  (45)  indicates  that 
2  “  )  “  P£(x,ct)  -  pQ  e~Z 


(46) 


and  hence,  since  z  <<  1  for  any  distance  of  interest,  the  wave  front  propagates 
is  such  a  manner  as  to  be  essentially  unattenuated  with  distance.  This  is 
in  contrast  to  the  case  where  forward  scattering  dominates  and  p^  drops  sharply 
to  zero  for  x  >  xc-  The  reason  for  this  difference  is  that:  a  spontaneously 
emitted  photon  which  is  created  near  the  wave  front  and  travels  in  the  forward 
direction  will  maintain  a  constant  distance  with  respect  to  the  front  and  hence 
upon  being  amplified  -  will  exhaust  the  laser  photon  intensity  at  that  point 
relative  to  the  wave  front;  whereas,  a  similarly  emitted  photon  which  travels 
in  the  backward  direction  will  rapidly  leave  the  region  of  the  wave  front  and 
hence  cannot  deplete  the  front  by  stimulating  emission.  Defining  a  quantity 


t 

o 


o 


(47) 


we  find  that  for  a  time  s  -  the  time  elapsed  subsequent  to  the  wave  front 
passing  a  given  point  x  -  for  s  >  t 


(u+l)e  -  y 


1  +  yz 


(48) 


That  is,  t  is  the  time  required  for  R  to  come  to  a  steady  state. 

O  X, 

Using  x  ■  i.  and  taking  Iq  to  be  the  maximum  intensity  which  can  be 

transmitted  through  the  atmosphere  under  the  assumption  of  negligible  back- 

2 

scatter  and  n  ■  1  (i.e.,  I  ■  124  kw/cm  ),  we  find  that  for  s  >  t 
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I  (£)  -*-4  kw/cm 

X 


For  this  same  choice  of  parameters 

t  ”  5  x  10  ^  sec 

o 

Hence,  the  steady  state  solutions  of  the  preceding  section  are  achieved  In  a 

time  of  the  order  of  t  •  5  x  10~^  sec.  If  we  wish  to  circumvent  the  various 

o 

limitations  and  difficulties  imposed  on  us  by  backscatter,  we  must  restrict 
ourselves  to  the  transmission  of  pulses  of  laser  light  with  durations  less 
than  5  x  10  ^  sec. 

It  is  of  interest  to  examine  the  interaction  of  a  laser  beam  with  a 
wave  packet  of  photons  at  the  first  Stokes  frequency  which  is  traveling  in 
tha  beam  but  antiparallel  to  it.  For  this  problem  it  is  most  convenient  to 
solve  the  initial  value  problem  , 

Pa(x,o)  -  L(x)  ,  pb(x,o)  -  B (x)  .  (49) 

Using  Eq.  (49)  in  Eqs.  (37)  and  (38),  we  readily  obtain 


(50) 


(51) 
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where 


r(x,y) 


2  (x_y)  +  2  \ 


dC[L(c5)  +  B(c0  3  • 


Let  us  examine  the  results  of  a  wave  packet  of  Stokes  waves  of  spatial 
width  W  and  photon  density  p^  colliding  with  a  laser  beam  of  photon  density 


L(x)  - 


B(x) 


That  is, 

let 

f 

for 

x  <_ 

0 

lo 

for 

X  > 

0 

/P1 

for 

o  < 

x  <  W 

lo 

Otherwise 

(52) 


The  -  wave  which  Is  generated  by  B(x)  lies  in  the  domain  o  <  t  <  W  in  the 
s  -  t  plane.  The  region  t  <  o  is  completely  unaffected  by  B(x) .  If  we  define 
the  gain  and  initial  value  parameters 


1  a. 

8i  "  2  +  afV  *  yi  "  ?  pi  » 


the  photon  densities  for  t  <  o  are 

1  + 


PjU.t) 


o  2g  t 
°o  , 
e  +  p 


(53) 


2g  t 


-  1 


pb(x*t}  ’  Po  ~2g  t 
o 

e  +  y 


(54) 
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Therefore,  in  the  absence  of  an  initiating  pulse:  the  back  scattered  wave  is 
tnonotonically  increasing  in  time  and  pfc  approaches  the  value  Pq  as  t  -  The 
laser  beam  intensity  is  monotonically  decreasing  and  eventually  decays  to  zero 

The  pulse  shape  is  described  by 


R^  (s  ,t) 


Pj^Vj+1)  +  P 

v  +  1  -  y_ 
o  o 


(55) 


The  wave  front  of  the  pulse  is  given  by  t  -  o.  At  a  time  which  is  long  compared 
to  (2go)_1  in  yQ  the  pulse  height  approaches 


^(-.o)  -  P1^0+1^  +  p0  "  p0(yl+1)  +  P1  ‘ 


(56) 


The  laser  beam  is  therefore  capable  of  amplifying  the  photon  density  to  a  level 
which  is  far  in  excess  of  Pq.  This  amplification  however,  is  limited  by  the 
spontaneous  depletion  of  the  laser  beam  dictated  by  Eq.  (53).  If  we  wish  to 
amplify  a  pulse  by  passing  it  through  a  cell  of  Raman-active  material  of 
length  d  we  will  generally  find  that  for  a  reasonable  set  of  assumptions  about 
the  parameters  available  in  the  laboratory 
d 


and  hence,  the  maximum  pulae  height  is  limited  by  the  length  of  the  cell  rather 
than  the  spontaneoua  creation  rate  B.  In  that  case,  the  height  of  the  emerging 
pulse  is  given  by 

d 

I 

V2  f  ■ 0>  *  °i e 


(58) 


The  half  width  of  this  pulse  is  given  by 


T 

O 


-op 


o  c 


e 


(59) 


which  will  generally  be  an  extremely  short  time.  The  ratio  of  the  emerging 
pulse  energy  to  the  energy  stored  in  that  part  of  the  beam  which  is  present 
in  the  cell  at  t  -  o  is  approximately 


U 


W 

d 


2c 

apod 


-ap 


W 
1  c 


e 


(60) 


On  the  basis  of  Eqs.  (58),  (59),  and  (60),  it  would  seem  to  be  feasible  to 

concentrate  virtually  the  entire  energy  in  a  laser  beam  into  a  pulse  with  an 

13  2  -12 

Intensify  of  the  order  of  10  watts/cm  and  width  of  the  order  of  10  sec. 

Experiments  have  been  performed^  in  which  the  initiating  pulse  appeared  to 

be  generated  by  Townes  Self  Focusing  at  the  end  of  the  cell.  The  emerging 

pulse  had  U  *  1/5  and  Tfl  a  3  x  10  ^  sec. 
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VII.  CONCLUSIONS 

We  have  found  that  Raman  backscattering  is  inconsequential  as  a  loss 
mechanism  in  the  propagation  of  a  beam  of  laser  light  vertically  through  the 
atmosphere  if  the  laser  light  is  composed  of  pulses  of  duration  no  greater  than 
5  x  10-5  sec,  or  if  a(o)  exceeds  a  (it)  by  eight  percent  or  more.  If  neither 
of  these  requirements  are  met  then  the  source  intensity  necessary  to  transmit 
the  maximum  amount  of  energy  through  the  atmosphere  would  be  increased  by  a 
factor  of  two.  In  that  case,  the  intensity  of  the  back  scattered  wave  at  the 
source  will  be  one-half  the  source  intensity  and  considerable  damage  to  the 
source  and  its  environs  may  result.  It  is  therefore  recommended  that  experi¬ 
ments  to  determine  the  asymmetry  of  the  Raman  cross  section  be  performed  if 
pulse  durations  in  excess  of  5  x  10  sec  are  desired. 

The  maximum  intensity  of  laser  light  which  can  be  transmitted  vertically 
through  the  atmosphere  without  frequency  change  due  to  Raman  scattering  is 
critically  sensitive  to  the  laser  line  width  (Av^/c) .  If  (Av£/c) 

1.2  x  10"3  cm-1,  and  if  back  scatter  is  inconsequential  then  for  a  beam  with  a 
10  cm  diameter,  this  maximum  intensity  is  27  kw/cm  ,  whereas,  if  Av£/c  is 
of  the  order  of  the  pressure  broadened  line  width  (9  x  10  cm  )  the  maximum 
intensity  is  120  kw/cm2.  It  was  suggested  that  intensities  far  in  excess  of 
this  figure  may  be  transmitted  in  the  form  of  a  higher  order  Stokes  wave.  If 
we  wish  to  transmit  energy  through  a  window  in  the  atmospheric  transmissivity 
of  width  Nv  (v  being  the  frequency  shift  due  to  Raman  scattering) ,  then  the 
maximum  emergent  N-th  order  Stokes  wave  will  be  approximately  N  times  the 
figure  cited  above.  However,  this  suggestion  did  not  take  the  probably  deleterious 
effects  of  defocusing  on  each  successively  higher  order  Stokes  wave  into  account. 

A  theoretical  investigation  of  this  three-dimensional  effect  as  well  as  a  more 
precise  evaluation  of  the  parameter  should  be  performed. 

Finally,  the  technique  for  producing  high  intensity  pulses  of  extremely 
short  duration  discussed  at  the  end  of  Section  VI,  appears  to  be  potentially 
quite  promising.  Further  investigations  along  these  lines  may  prove  to  be 

rather  fruitful. 
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